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Abstract 

This article proposes the construction of Wigner measures in the infinite dimensional bosonic 
quantum field theory, with applications to the derivation of the mean field dynamics. Once these 
asymptotic objects are well defined, it is shown how they can be used to make connections between 
different kinds of results or to prove new ones. 
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1 Introduction 

The bosonic quantum field theory relies on two different bases : On one side the quantization of a 
symplectic space, the approach followed for example by Berezin in HBerL Kree-Raczka in [KrRa); on 
the other side the gaussian stochastic processes presentation also known as the integral functional point 
of view followed for example by Glimm-Jaffe in MGlJal and Simon in HSimll . Both approaches have 
to be handled in order to tackle on the most basic problems in constructive quantum field theory (see 
[BSZ][DeGe]). The interaction of constructive quantum field theory with other fields of mathematics 
like pseudodifferential calculus (see [BeShl or MLasID or stochastic processes (see |Mey[ [AtPa]) is often 
instructive. 

In the recent years the mean field limit of A-body quantum dynamics has been reconsidered by various 
authors via a BBGKY-hierarchy approach (see HESY11 IIESY21 llFGSl llFKPll IBGGMI |Spo) and ISerl for 
a short presentation) mainly motivated by the study of Bose-Einstein condensates (see HCasll "). Although 
this was present in earlier works around the so-called Hepp method (see [HepJ and |GiVe]), the relation- 
ship with the microlocal or semiclassical analysis in infinite dimension has been neglected. Difficulties 
are known in this direction : 1) The gap between the inductive and projective construction of quan- 
tized observable in infinite dimension; 2) the difficulties to built algebras of pseudodifferential operators 
which contain the usual hamiltonians and preserve some properties of the finite dimensional calculus 
like a Calderon-Vaillancourt theorem, a good notion of ellipticity or the asymptotic positivity with a 
Garding inequality; 3) even when step 2) is possible, no satisfactory Egorov theorem is available. 
Recall the example of an A-body Schrodinger hamiltonian 

H N = -A+^ £ V( Xi -Xj) , onR dN , 

l<i<j<N 

and consider the time-evolved wave function 

y N (t) = e - itHN y® N , ¥ e L 2 (R d ) . 

The 1 -particle marginal state, the quantum analogous of the one particle empirical distribution in the 
classical A-body problem, is given by 




The mean field limit says that in the limit N — > °°, the marginal state evolves according to a non-linear 
Hartree equation 

p\t) = \z(t))(z(t)\+o(l) , astf-x», 

with { id tZ = - A z+( v *\z\ 2 )z onR t xW l 
\ z(t = 0) = y. 
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By setting N = ^ and in the Fock space framework with £ -dependent CCR (i.e: [a(g),a* (/)] = e (g, /}), 
the problem becomes 

1 



1 

Hn — - 
s 



Va*(x)Va(x) dx+ I V (x — y)a* (x)a* (y)a(x)a(y) dxdy 

'■ Jm 2 * 



s 



e ~itH N _ e -ijH e 



Tr [Ap\t)] = (W N (t) , dT(A)¥ N (t)) = (W N (t) , p A (z) Wick W N (t)) , 

where p A is the polynomial Pa(z) = (z,Az) . Higher order marginals, taking into accounts correlations, 
can be defined after using the polynomials p A (z) = {z® k ,Az® k ) with A e ^f(L 2 (R kd )) . 
On this example, the scaling of the hamiltonian, of the time scale and of the observables as Wick 
operators enters formally in the £-dependent semiclassical analysis. The Hepp method concerns the 
evolution of squeezed coherent states ( [Hep 1 1 Gi Vel MCasIl ). which amounts in the finite dimensional case 
to the phase-space evolution of a gaussian state according to the time dependent quadratic approximation 
of the non linear hamiltonian, centered on the solution to the classical hamiltonian equation. We refer 
the reader to [CRR] for accurate developments of such an approach in the finite dimensional case. 
In the nineties and as a byproduct of the development of microlocal analysis, alternative and more 
flexible methods were introduced in order to study the semiclassical limit with the help of Wigner (or 
semiclassical) measures (see [Bur][Ger][HMR][LiPa][Tar]). Such objects are defined by duality and 
rely on the asymptotic positivity of the £-dependent quantizations. It gives a weak but more flexible 
form of the principal term of the semiclassical (here mean-field) approximation. Via the introduction 
of probability measures on the symplectic phase-space, it provides an interesting way to analyze the 
relationship between the two basic approaches to quantum field theory. Further in finite dimension, the 
Wick, anti-Wick and Weyl quantizations are asymptotically equivalent in the limit e —> 0. This is not so 
obvious in infinite dimension. 

Several attempts have been tried to develop an infinite dimensional Weyl pseudodifferential calculus 
with an inductive approach. Lascar in [Las] introduced an algebra and a notion of ellipticity in this 
direction, making more effective the general presentation of [KrRa]. The works of Helffer-Sjostrand 
in jHel2][HeSj] and Amour- Kerdelhue-Nourrigat in [AKN] about the pseudodifferential calculus in 



large dimension motivated by the analysis of the thermodynamical limit enter in this category. With 
such an approach, it is not clear that the infinite dimensional phase-space is well explored and that 
no information is lost in the limit £ — > 0. Meanwhile this inductive approach is limited by Hilbert- 
Schmidt type restriction like in Shale's theorem about the quasi-equivalence of gaussian measures. It is 
known after IIGrol that the nonlinear transformations which preserve the quasi-equivalence with a given 
gaussian measure within the Schrodinger representation are very restricted and do not cover realistic 
models. Hence no Egorov theorem can be expected with Weyl observables. 

Simple remarks suggests alternative point of views. The Wick calculus with polynomial symbols present 
encouraging specificities: It contains the standard hamiltonians, it makes an algebra under more general 
assumptions (the Hilbert-Schmidt condition can be relaxed) and allows some propagation results when 
tested on appropriate states (see [FGS][FKP]). Meanwhile the Wigner measures in the limit £ — > can 
be defined very easily via the separation of variables as weak distribution, in a projective way which fits 
with the stochastic processes point of view. 

After reviewing and sometimes simplifying or improving known results and techniques about the mean 
field limit, our aim is to show the interests of the extension to the infinite dimensional case of Wigner 
measures: 

• After the introduction of the small parameter e —* and the definition of Weyl operator W(z), 
zG 5° the phase-space, choosing between the quantization of symplectic space and the stochastic 
processes point of view is no more a question of general principles nor of mathematical taste. 
It is a matter of scaling. The symplectic geometry arises when considering macroscopic phase- 
space translation W(|), while the operator W(z) is used with this scaling in the introduction of 



3 



Wigner measures via their characteristic function. Corrections to the mean field limit considered 
for example in MCCDB with a stochastic processes point of view can be interpreted within this 
picture: They attempt to give a better information on the shape of the state in a small phase-space 
scale. 

• Once the Wigner measures are well defined as Radon measures, it is possible to make explicit the 
relationship between different kinds of results and to extend them in a flexible way. It accounts 
for the propagation of chaos (result obtained via the BBGKY approach) according to the classical 
hamiltonian dynamics in the phase-space. Actually we shall prove in a very general framework 
that the propagation of squeezed coherent states as derived via the Hepp method implies a weak 
version of the mean field limit for product states. Further propagation results can be obtained for 
some non standard mixed states without reconsidering a rather heavy analysis process. 

• The comparison between the Wick, Weyl and anti-Wick quantization can be analyzed accurately 
in the infinite dimensional case. With the Wick calculus, complete asymptotic expansions can be 
proved after testing with some specific states. The relationship of such results with the propagation 
of Wigner measures works in a rather general setting but has to be handled with care. 

• The gap between the projective and inductive approaches can be formulated accurately in the limit 
£ — ► 0. We shall explain in the examples the possibility of a dimensional defect of compactness. 

This work is presented and illustrated with examples simpler than more realistic models considered in 
other works like MGiVe 1 [Hep 1 [ ES Yll llESY2l HBGGM1 with more singular interaction potentials. That 
was our choice in order to make the correspondence between various approaches more straightforward 
and to pave the way for further improvements. We hope that this information will be valuable for other 
colleagues and useful for further developments. 

The outline of this articles is the following. In Section[2l standard notions about the symmetric Fock 
space are recalled and Wick calculus is specified. In Section [3] the Weyl and Anti-Wick calculus are 
introduced in a projective way after recalling accurately (most of all the scaling) of finite dimensional 
semiclassical calculus. The Section @]recalls the distinction between coherent states and product or Her- 
mite states, and their properties when measured with different kinds of observables. The two methods 
used to derive the mean field dynamics, the Hepp method and the analysis through truncated Dyson 
expansions, are reviewed within our formalism and with some variations in Section [5] The Wigner 
measures are introduced in Section [6] with the extension of some finite dimensional properties and spe- 
cific infinite dimensional phenomena. Finally examples and applications are detailed in Section |7J in 
particular: 1) reconsidering a simple presentation of the Bose-Einstein condensation shows an interest- 
ing example of what we call the dimensional defect of compactness; 2) a general result says that the 
propagation of squeezed coherent states, which can be attacked via the Hepp method, implies a slightly 
weaker form of the propagation of chaos (formulated with product states and Wick observables); 3) the 
mean field dynamics can be easily derived for some states which present some asymptotically vanishing 
correlations. 

Acknowledgements: The authors would like to thank V. Bach, Y. Coudene, J. Frohlich, V. Georgescu, 
C. Gerard, P. Gerard, S. Graffi, T. Jecko, S. Keraani and A. Pizzo for profitable discussions related with 
this work. This was partly completed while the first author had a sabbatical semester in CNRS in spring 
2007. 

2 Fock space and Wick quantization 

After introducing the symmetric Fock space with £ -dependent CCR's, an algebra of observables result- 
ing from the Wick quantization process is presented. 
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2.1 Fock space 

Consider a separable Hilbert space 3: endowed with a scalar product (.,.) which is anti-linear in the 
left argument and linear in the right one and with the associated norm \z\ = \J (z,z). Let a = Im(., .) 
and 5 = Re(., .) respectively denote the canonical symplectic and the real scalar product over 3f. The 
symmetric Fock space on 3? is the Hilbert space 

jf = ^\j3f = T s {3f), 

n=Q 

where V" ^ is the n-fold symmetric tensor product. Almost all the direct sums and tensor products are 
completed within the Hilbert framework. This is omitted in the notation. On the contrary, a specific alg 
superscript will be used for the algebraic direct sums or tensor products. 

For any n € N, the orthogonal projection of % n 3? onto the closed subspace V" will be denoted 
by y n . For any (&,&,...,&,)€ 2™, the vector & V £ 2 V • • • V $ n G V" & will be 

§lV§ 2 V"-V^ l = ^ 1 (§i®&"-®§,) = ^ £ &7(l)®§ff(2)-"®&j(n) 

The family of vectors (£i V • • • V £,n)fe%- is a generating family of \J n ' al ^ 3f and a total family of V" 3f . 
Thanks to the polarization identity 

L n - e,=±l 7=1 

the same property holds for the family {z® n ) ne ^ ze ^ ■ 

For two operators A k : \J ik 3f -> V A ^, A: = 1,2, the notation A\ \J A 2 stands for 

h+h ji+h 

A l \jA 2 = y h+h o{A l ®A 2 )oy h+h e^{\J 3f,\J 3f). 

Any z£ 3f is identified with the operator from \J Q 3f = C 3 X ^ Xz £ 3f = \J l 3f while (z\ denotes the 
linear form 3f 3 ^ i— > (z, £;) G C. The creation and annihilation operators a*(z) and a(z), parameterized 
by £ > 0, are then defined by : 

a {z)\\j n s = Ven {z\®Iyn-\ 

a*(z)\y n 2° = V e ( n +1) ^n+i°{z®I\i n &) = ^Je{n+\)z\Jl\i"^. 

Each of (a(z)) ze j° and (a*(z)) zG ^ are commuting families of operators and they satisfy the canonical 
commutation relations (CCR): 

[a(z 1 ),a*(z 2 )} = e{zi,z 2 )I. (2) 

We also consider the canonical quantization of the real variables <I>(z) = -^( a *(z) +a(z)) and H(z) = 
&(iz) = ~j^{ a {z) — a *(z))- They are self-adjoint operators on J4? and satisfy the identities: 

Mz l )Mz2)] = iea(z l ,z 2 )I, [<fr(zi),n(z 2 )] = ieS{zi,zi)L 
The representation of the Weyl commutation relations in the Fock space 

W(zi)W(z 2 ) = e-'i a ^'^W{zi+z 2 ) (3) 
= e- iea ^^W(z 2 )W(zi), 
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is obtained by setting W (z.) = e'®^ . The generating functional associated with this representation is 
given by 
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(fl,W(z)Q>=«- 

where £2 is the vacuum vector (1,0, •■■) G Jff. The total family of vectors E(z) = W f-^rj ^ 
e \[a*{z)-a{z)\Q z e have the eX pii c i t f orm 



E(z) 



_Hi f 1 a*(z)" 
n=0 fc n ■ 



e 2e 



-n/2_ 



(4) 



«=0 vm 

The number operator is also parametrized by e > 0, 

N\ y»sr = £n/|y» jr. 
It is convenient to introduce the subspace 

alg n 

neN 

of J^, which is a set of analytic vectors for N. 

For any contraction S G J*?(3f ), \S\^r^ < 1, r(5) is the contraction in defined by 

r(5)| V »^=5(8)5---(8)5. 

More generally can be defined by the same formula as an operator on Jtfi„ for any B G JSf (JT). 
Meanwhile, for any self-adjoint operator A : ,3° D &{A) — ► 2, the operator dT(A) is the self-adjoint 
operator given by 



k 



For example = dT(I) . 



2.2 Wick operators 

In this subsection we consider the Wick symbolic calculus on (homogenous) polynomials. We will show 
some product and commutation formulas useful later for the application. For example time evolved 
Wick observables can be expressed as £-asymptotic expansion of quantized Wick symbols. For a de- 
tailed exposition on more general Wick polynomials we refer the reader to IDeGel. 

A {p,q) -homogeneous polynomial function of z G ^ is defined as Pe(z) = (-{z® q ,z® p ), where I is 
a sesquilinear form on (® qMg 2) x (<g> p > a,g 2?), with P t {Xz) = l q X p Pt{z). Owing to the polarization 
formula (Q} and the identity 

Jt) Jo 
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the correspondence £ *—>■ Pg is a bijection when the set of forms is restricted to the sesquilinear forms 
on (\J q ' alg 3f) x (\jP' alg 3f). Any of the continuity properties of Pg are thus encoded by the continuity 
properties of the sesquilinear form I with the following hierarchy (from the weakest to the strongest) 

|£(ij 1 V...Vtfe,§iV...V§ 1 ,)| <Q\r} 1 \ 3r ...\r] q \^\^\ 3r ...\^ p \^, tyG JF^-GiF 

\ltt,V)\<Ct\$\vgty\y*z> W^\J eV ^ © 

, *T g N, Cy g C, 



£ Cy^-,V^)I<Q 

i<U<K 



(V f ^)*®(V p 50 



For example, when = g = 1 the two first ones define J£(3f), while the third one defines the space of 
Hilbert-Schmidt operators. By Taylor expansion any (p,q) -homogenous polynomial P admits Gateaux 
differentials and we set 



d£d*P(z)[ui,~- ,u k ,vi,-~ ,vk] = d ui ■ ■ ■ d Uk d n ■ ■ ■ d Vk ,P(z) 
where d u ,d v are the complex directional derivatives relative to u,v G 3f. 

Definition 2.1 For p,q G N, the set of (p ^-homogeneous polynomial functions on 3f which satisfy 
the continuity condition (0) is denoted by £P pq (3f): 

[b(z) G ^ M (^)j ^ | ^ = 

77ze subspace of & p q(3f) made of polynomials b such that b is a compact operator b G Jf°°(J2?) (resp. 
b G ST(3T)) is denoted by ^ q (3f) (resp. @> r p q (2?)). 

It will be sometimes convenient to consider b as an operator from P 3f into 0* 3? with the obvious 
convention for symmetric operators b = S? q bS? p . Owing to the condition b G Jz? (V p -3°, V 9 f° r 
b G ^ p>q (3f), this definition implies that any differential d^d^b(z) at the point z G 3f equals 



(6) 



We will mainly work with fixed homogeneity degrees but the key statement of this section (Propo- 
sition |2]6l) says that ©p^eN'^/'^-^') is an algebra of symbols with the same explicit product formula as 
in the finite dimensional case. 

With any "symbol" b G & P:q (3?), a Wick monomial b w,ck can be associated according to: 



with£ = (p!)- 1 ^!)- 1 ^^^) . 

Here are the basic symbol-operator correspondence: 

<— >/2S(£,z) <— <&(§) (z,Az) <— rfT(A) 

(^z) — a(§) V2a(^,z) <— > n(£) |z| 2 <— > AT. 

Other examples can be derived from the next propositions. The first one is a direct consequence of the 
definition f7]). 
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Proposition 2.2 The following identities hold true on M'finfor every b G & p q {3f): 

(i) (b wick y =h wick . 

(ii) (C(z)b(z)A(z)) Wick = c Wkk b Wick A Wick , if A G 0> a ^S), C G 5»o,p(^)- 
e I 'i dr ( A )fo w, ' cfe e -''l dr ( A ) = (i,(r M z)) WH if A is a self-adjoint operator on 3f. 

p i 

Proposition 2.3 (i) The Wick operator associated with b(z) = Y[ (z, X Y\ (^j^)> Hi' %j ^ =^°> eaua ^ s 

i=l j=\ 

b Wick = a*(T ]l )-a*(Tlp)<Zi)-a(Z q y 
(ii) For b G & > p^(3f) ant/ z.& 3f the equality 



(z® i ,b Wick z® k ) = 8+ pJ _ q 



k\j\ 



(k-p)\{j-q)\ 



\z\ k ~ p+hq 



b{z) 



(8) 



holds for any k,j G N. 77ie symbol 8^ denotes 8 a pl^ + ^(a) where 8 a p is the standard Kronecker 
symbol. 



□ 



Proof, (i) is a direct consequence of Proposition |Z2] with ((z, £)) Wick = a* {£,) &n& ({£, ,z)) Wick = a(E,) 
(ii) This comes directly from the definition (0 of b Wlck . 

The next result specifies the boundedness properties of b w . 
Lemma 2.4 For b G & PA (3f), the estimate 



b w \^^^^)^ 5 t P ,j- q ue) 2 m^\b\ Nq with t=—dpib, 

y -H- 



holds for any k,j€N. 
This implies 



(N)—2b Wick {Ny? 



< 



(9) 



(10) 



Proof. A consequence of ® is b Wick ( V k 3T) C V ; 3T with j = k- p + q. For y G V & and j = k- p + q, 



write 



\b Wick vL 



i\k\ 



V - (k-p)\ 

< (J£) f (te)^ 



yjib^I^k-p^Xjr 



k\ 



(j-qV-P V (k-p)\kP 



IVlv*jr 



□ 

An important property of our class of Wick polynomials is that a composition of bY o fe!p c * with 
&i,&2 G ®pfq&p,q{ 3f) is a Wick polynomial with symbol in ®p%& P ,q{ 3f). In the following we prove 
this result and specifies the Wick symbol of the product. 
For b G &p^ q (3f), specific cases with j = or k = of © imply 



d k b{z) G (\/ 3f)* and dlb(z)e\j2f, 

for any fixed z G 3f. For two symbols G i = 1 , 2, and any k G N, the new symbol di^i .d|&2 

is now defined by 



d k b 1 .d!b 2 (z) = (d k b l (z),d!b 2 (z)) 



(V* a?)* ,v* % 



(ii) 
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We also use the following notation for multiple Poisson brackets: 

{b u b 2 }^ = dzh-dlh- d k z b 2 .dlb x , 
{biM = {hM m - 

These operations with polynomials are easier to handle than there corresponding versions for the op- 
erators hi G .if (V p ' V* ^H- Nevertheless their explicit operator expressions as contracted products 
allow to check that ©p^^p ,q{^) is stable w.r.t these operations . 

Lemma 2.5 Fix pi,p 2 ,q\ and q 2 in N. For two polynomials b[ G 3P piAi {2£), i = 1,2, set hi = (p,!^,!) -1 
df'dfbi and for any k G {0, . . . ,min{/?i ,172}} 



h Qh 2 = — — J-— — -af +*-*af + * 2 ~* 

(Pi +P2-^)!(<?i +<?2-^)! 



Then 



pi+p 2 -k q l +q 2 -k 



with the estimate 



_ it _ 
&i©&2 



— / 7 \ 1 / im l v(\ip\ <¥ \m <v\ "2 a>(\/P2 v \m >¥\ ■ v^-V 



Proof. For y G V Pl ^° and G V ?2 ^> introduce the vector 

'^d*b d (z)e 1 

with &0 (z) = (z 92 , 0) and the form 



<z»*-*,*> = ((z^- k \ 8/3*,) = (®-*)!# M z) G v^r 



(^^-^^^-^^(^G^irr, with Mz) = (iff,^). 

The identity 

((^z^-^),^^^))^^^^ = (t^^z^-^^^-^^)^^-^ (14) 

is obviously true when i/a = ^®P l and = r\® q2 with ^,T] £ 3f. Since (£,® n )E e % is a total space of 
V" ^ with the polarization identity ©, the identity d) holds for all G V 92 & and all i//" G V Pl ^ ■ 
After noticing the relations 

with y/- = ^z® 91 and <j> = h 2 z® Pl , the identity £[4]> leads to 

Therefore d^b\.d\b 2 is a continuous homogenous polynomial in 3 g p l +p 2 -k,q l +q 2 -k{S ? ) with the associ- 
ated operator given by (fT2l ). The estimate (fT3T > follows immediately by (fT2l ). □ 
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Proposition 2.6 The formulas 



(!) 



■jWick foWick 



'minjp!,^} £ k 



Wick 



k=0 



^^1(2)62(0)1^) 



ll'/c* 



^max{min{pi,^ 2 },miii{p2,?l}} g jt 



Wick 



{bub 2 } {k) 



(15) 
(16) 



hold for any bj £ B^p^q.^Sf), i = 1,2 as identities on J^fi, 



Remark 2.7 T/z/s resw/f /jas exactly the form of the finite dimensional formula. Lemma 12.51 gives 
relation with the writing which can be found in HFKPV . 

Proof. The second statement (ii) is a straightforward consequence of the first one (i). Let us focus on 
(i) which will be proved in several steps. 

Step 0: Before proving the identity, first notice that both sides are well defined. Actually, for any 
b £ £Pp !q (3f), the operator b Wlck sends M'fin into itself. Hence, the product b^ ,ck oti% lck is well defined 
as an operator — > Finally we know from Lemma [231 that e e ^ z,dm 'b\ (z)b2(co) | _ ffl belongs to 

m al g <«l / <*"\ 
*X>p,q<^ r p,q\-^- )■ 

Step 1: Consider &i(z) = (77 , z) and 6 2 (z) = (2, € N. The formula 



fl (T7K(?) 9 = fl*(§)Vn) + e?<i?^>fl*(« 



W- 1 



is exactly 

b mck b mck = {p ib2 fick + e ( dzbl _ d . bl) 

Step 2: Consider 61(2) = P P (z) = (ij , z) 7 * and 62(2) = (2, G N. The induction is already 

initialized for p = 1 according to Step 1. Assume that the formula is true for p — 1 and all q G N and 
compute 



WickuWick 



pmck b 



pWick pWickjjWick _ ^W 



mm{p-l,q} £ k 



k=0 



Wick 



fl(Tj) 



min{p-l,<?} it 



e\ ~k_q\ (P-1)1 ^ ( g )g -* fl(l?r i-t 



it! 



min{p-l,g} & 

fcP (o-jfc)!(/?-l-jfc)! 



I 

>t=0 



' ( 9 -*)! (p-1-*)! 

q\(p-\)\ 



a*(^- k a(v) p - k 

+e(q-k)(l 1 ^)a*^- k a(riy- {k+l) 



e k ^^Yg\(p-\)\ 
- *!( 9 -*)!(p-l-*)! 



a*(^)«- fe a(r])^* 



min{p,q} Q k 



'd%,d k b 2 



k=Q 



kl 



Wick 



We used several times the relation 



n\ 



(n-j)\ 



(■n,z) n - J (n\ 9J 
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and its dual version for d[b2 ■ 

Step 3: From Step 2, the statement (ii) of Proposition [2]2] leads to 



a*(Z 1 ) qi a(ri 1 ) Pl a* (<§ 2 ) 92 a{r} 2 ) P2 



min{pi,g 2 } c k 



k=0 



kl 



Wick 



for any ^ 1 ,^ 2 ,T] 1 ,T] 2 S 3f and any p\,q\,P2,qi £ N. Again the polarization formula (Q]) in the form 



1=1 



yields the result for any 



that is for any in the form 



pi 



e,=±l 



7=1 



M«)=n<*.#>IIft;>*>» £ = 1 ' 2 ' 

7=1 



i=\ 



^ = |^V...vO(T?fv...V7]l| 5 £=1,2. 



(17) 



Step 4: We want to check the identity 



Wn' ,1>1 °b 2 Wn 



£ -U n ,,(d z p hdib 2 ) mck ¥n 



p=0 



pi 



for any fteV"^° and any y„' G V" ^,n,«'el 
From the definition of b Wlck , the left-hand side equals 

Wn' > bi lCk ob2 Wk Xjf n \ = C„,„' jPl2 , 9l2; e ^VA)' , (bl\J l\y„+ n - P2 -p l ^) (bl \ll\y- P1 s ) % 

= C n>n ' )Pl2iqi2)E ((t>i V^|y'-«i ir) VV > (^V^l v"~''i 3") 
Similarly and owing to Lemma [231 every term of the right-hand side satisfies 

^,(dibi^b2) Wkk Yn 



: C'n.n 1 ,p,pi2,q\.2,e (^"' ' (^1 ®-*®»-*.yJ ( 



> liwrx-p <sjf} [ IffiPi-P J" 8> £>2 J \J I\/"-p\-P2+P $f 



Wn 



r' 



n,n' ,p,p\2 ,<?i.2,£ 



1 ®^2)"'-''i SJ V n ' ' (/(g>'' 1 ~''3* ® ^ 2 ® I\j"-Pi-Pi + P ) Wn) ■ 



Hence for fixed WniWn 1 G ^/m> both side are sesquilinear continuous expression of (£1,62) when the first 
factor is considered with the *-strong topology of operators and the second one with the strong topology. 
The operators (fTTT ) for which the equality is true, form a total family for these topologies: In two steps, 
approximate first any finite rank operators and then bounded operators by finite rank operators. Thus 
the equality holds for any b e G ^ PiM (2?), £ = 1,2. □ 

Remark 2.8 The formulas (1151 ) and (1161 ) make sense with e-dependent symbols. One can work with 
polynomials in £ 

b( z ,e)=j^e a b a (z), b a e<? M (2?) 



a=0 
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or with asymptotic sums 



b{z,e)~ £e a fc«(z) ^eV). 



a=0 



/7ie expression (1151 ) a«<i (1161 ) ta/ce ?/ze?i the form 



Wick 



jWick^Wick 



iWick ,Wick 



7=0 \a+/3+*:=; 



j=l \a+p+k=j 



for by ~L a £ a »i ; a G ^piaii^O and b 2 ~ Y,p S p b 2 fi G ^^(JT) . //ere (rasp. (p 2 ,q 2 ))does 

not depend on a ( resp. j8 J. 

We have the following useful result. 

Proposition 2.9 For any b G ©^n^ 2 /^^) w<? have: 

( i) b Wwk is closable with 

M = vect{W(z)0,0 G J^ fin ,z G iT} 

a core of the closure. 

(ii) By setting E{z) = W(^r)£l according to ©, the identity 



b(z) = (E(z),b mck E(z 



holds for every z G £tf . 

( Hi) For any zq G 3? the identity 



(18) 



W(^zoTb Wick W(^zo) = (Hz+zo)) wick 

IE IS 
holds on J% where b(-+z ) G ®% en ^ P , q {^) ■ 

Proof, (i) b Wlck is closable by Proposition [272](i). It is enough to consider b G & p , q (3?) when we prove 
that is a core for the closure of b w,ck . The last statement is deduced from the estimate 



1 



>(*)! 



£ b wick^ {zr(p (k) <\~by {yP ^ yqv \<p^\ yk ^x 
«=()"• M 



n=0 



kl 



for any <pw £\J k 3f and z G 3f. In order to prove (fl9l . use Lemma l2T4l and estimate the action of b Wwk on 

Wicki 



0>(z)"<pW by max |^^' c *|jg'(y and bound the norm of 0(z) n <pW by |<p w 1 

(ii) One writes for o G ^ P , q (T) and z G J° 



(*)|| 7 |n . / (2£)"(n+<:)! _ 



(/?(z),& w,v */?(z)) = e* £ 



fyWick^(g)n2\ 



/n\\\Jn 2 \ 



_ Elf 



£ 2 Z 



ni-p+n 2 -q 



»i ■/.':■-• ' ''' P -?) ! \/( n 2 -/>)!£ 
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, al.2 



J p,qeN^ PA 



& P: q(3f) come from the Taylor expansion and 



(iii) The fact that b(.+zo) remains in the class 

©. In order to prove the equality, differentiate A(t) = W(^§tz )b(z + tz ) Wick W(^§tzo)* 
sense on J%. Proposition 12. 61 implies 



in a weak 



id t A(t) 



,V2 x 
W — tzo) 

IS 



,V2 



- M^—Zo), Hz + tzo) mck ] + id t b(z + tzo 



.Wick 



W^tzoT 

IS 



wC§tzo) [(izo,d- z b(z + tzo)) - (d z b(z + zo) , iz )+id t b(z + tz )] Wick W(^tzoy = 0. 



□ 

Remark 2.10 The relation (1181 ) allows to define easily the Wick symbol of an operator which is defined 
as a series, when it makes sense, instead of a Wick polynomial. For example the Wick symbol of the Weyl 
operator W(%) equals 



(E(z),W(£)E(z)) = (a,e- iea £>^W(Z)Q.)=e iV * s ^) e -^ L . (20) 

A variation of Proposition |2.9| ensures that b(Az + Zo) can be Wick quantized for any bounded complex 
affine transformation in 3? when b E ^ Ptq (3f). Actually real symplectic affine transformations of sym- 
bols in & p ^ q (£t) may also be Wick quantized but only under a Hilbert-Schmidt condition on A which 
agrees with Shale's theorem or the presentation of general Bogoliubov transformations (see IIBerlO . The 
following result will be useful in Subsection 15. II 

Proposition 2.11 Let B E J£(3f) and let B 2 E J2? 2 (iF) be an Hilbert-Schmidt operator on 2 and let 
J : 3- 9 z 1— > Jz='Z E 3f be any anti-unitary operator on 2. Then for any b E ZP p>q (£>f) the polynomial 
b(Bz + B 2 z) belongs to ®pi+ q i =p + q &p\ q '(2) with the estimate 



d^dP'b(Bz + B 2 z) 



<Cp, q {\B\^ {S) + \B 2 \ 



\p+q . 



^(\/ p 3f,\Ji 3f) 



Proof. For b E Pp^SF) write, after recalling b = ^ q l^ p in «5f (<g) p 3T, <g)« St), 
b(Bz + B 2 z) = ((Bz+B 2 zf q ,b(Bz+B 2 zf p ) 

= L £ C i C P ( ( B z) 0y - j ® (Bizf j , b(Bizf k ® (Bz)® p - k 

j=0k=Q 

= tt C J q C k £j tk (z^ +k - j , z® p+j - k ) . 
j=ok=a 

The sesquilinear form l j;k is defined on ((g)^' 3? ® alg ®* 3f) x (® J 2 ® alg ® p ~ k 3T) by 



j,k (0i ® <h, ¥1 ® Y2) = ( (B®"-^! ) 03) (fif W) , b(Bf<h) ® (B 



®p—k^ 



Vl 



It satisfies for = £^ =1 <j) ha <8> <h, a and = £^ =1 ® Vifi 

£jt(*,V) = t({BfW^),C^ p - k ) Wlfi 
p=i 

= t(wpAB* 2 r j c^ p - k ) ¥up 

13=1 



N 
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with 

N P~k 1 

a=l \ / 

Since B^ J is a Hilbert-Schmidt operator the estimate 

holds for any *P G (g) 7 ® al z <g) p ~ k 3f . In order to estimate IC^I^,^-* ^. ^ 6J) take any U e 3f 
and any V G <8) p_ * J° and compute 



\(U,C*V)\ 



JT (h, a ,(B*)®*-JC uv B$ k h ta 



a=\ 



with Q/ V = (/®f-y^®<^|)*(V^®|V» e-§f((g)ir,(g) 
Again the Hilbert-Schmidt condition implies 

\(U ,CcpV)\ < \Bz\jf2^ \ B \%>(S) \ u \&2r m^wsr^S) \ V \&- k % '1*1 (S'-^jr ■ 
We have proved an estimate for \C&\ which implies that the estimate 

\^j,k(^^)\ - I-^I^jq l^lj^jT) ^ |^|jSf(Vf ^V 4 ^) l*l® 4_y+ * ^ ' 

extends continuously to any <I> G (g) 9- - 7 ' -1- * and any *F G Jf. It holds in particular when <& G 

yq-j+k g and y G Hence G J^) holds for any (y,Jfe), j < <? and * < p, 

with a norm estimate which yields the final result. □ 

3 Weyl and Anti-Wick quantization 

Our extension of the Weyl and Anti-Wick pseudodifferential calculus to the infinite dimensional case is 
based on a separation of variables approach within a projective setting. This is slightly different than the 
one developed by B. Lascar in ULasI where the inductive approach leads to a natural Hilbert-Schmidt 
condition and restricts the exploration of the infinite dimensional phase-space 3f. 

3.1 Cylindrical functions and Weyl quantization 

Let P denote the set of all finite rank orthogonal projections on 3; and for a given p G P let L p (dz) 
denote the Lebesgue measure on the finite dimensional subspace p3f. A function / : 3f — > C is said 
cylindrical if there exists p G P and a function g on p3f such that f(z) = g(pz), for all z G 3f. In this 
case we say that / is based on the subspace p3f. We set y cy i(3f) to be the cylindrical Schwartz space: 

(/ G y cy l{2?)) & (3p G P, 3g G y( P 3f), f{z) = g(pz)) . 
It is well known that the Fourier- Wigner transform defined by the expression 

z ^Y[ ( j ) , v ]( z ) = ( v ,W(V2Kz)<l>), 



14 



for any (j), \ff G Jt, belongs to L 2 (p£tf,L p (dz)) C\Co(p3f) for every pgP. Introduce the Fourier trans- 
form of a function / G o^ cv /(iF) based on the subspace as 



Jp2? 



and its inverse Fourier transform is 



J p3f 



Therefore the so-called Wigner transform can be written asW[^,y/]=^ 1 y]]- With any symbol 
b G y cy i(£¥) based on pif, a Weyl observable can be associated according to 



• Weyl 



pSf 



&[b]{z) W(V2nz)LJdz). 



(21) 



It can be expressed as a quadratic form in the following way 



(Y,b Weyl $) 



3te 



p2? 



&[b](z) n<l>M(z)L P (dz) 
b(z) W[$M(z)L p {dz). 

pS 

Note that b^ 1 is a well defined bounded operator on for all b G y cy [(3f) since Y[^,\jf](z) is a 
bounded function and <&[b] (z) is in L 1 (p^ ,L p (dz))- Remember also that this quantization of cylindrical 
symbols depends on the parameter e like the Weyl operators W(V2nz) ■ 

The next estimate will be useful. A similar inequality can be found in [DeGe]. 

Lemma 3.1 For any 8 G [0, 1] there exists a constant Cg > such that the estimate 

[W(zi)-W(z 2 )](N+\)- s / 2 \ <C S \zi-z 2 \ S [mm(e\zi\,e\z 2 \) s + mzK(l,e)% 

holds for all £ > 0, and all zi,z 2 G 3? . 

Proof. We have by Weyl's relation 

[W(zi) -W(z 2 )](N+ l)- 5 / 2 | < \[W(zi -z 2 )-I](N+ l)- 5 / 2 | + l^fe*) - 1 

The estimate \e' s — 1| <C§\s\ s , leads to 



(22) 



g fe<T(zi,a)_ j = g ie<T(zi-?2,a) _ j 



,i'ecr(zi,Z2-zi) _ j 



<C 5 £ 5 |zi-z 2 | 5 min(|zi|,|z 2 |) 5 . 



The first part of the r.h.s. in (1221 is estimated via a complex interpolation argument. Indeed, for 5 = 
notice that \W(zi — z 2 ) — 1\ < 2 and for 5 = 1 the estimate |e" — l| < Ci|s| combined with the spectral 
theorem yields 



[w(z 1 -z 2 )-/](w+ir 1/ v 



< Ci 

< Ci 



I*(zi-z2)i(iv+i; 



-1/2 



< J > (zi-z 2 )(^+l)" 1/ V 



Now by the number estimate (PTOl ) we obtain 

[W(zi-z 2 )-/](A^+l)" 1/2 < C max(l,e) \z\-z 2 \ 



□ 
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3.2 Finite dimensional Weyl quantization 

The finite dimensional Weyl calculus provides us a collection of results on the Weyl quantization. We 
specify here the relation between the Weyl quantization defined on Ht via (|2TT) and the usual semiclas- 
sical Weyl quantization within the Schrodinger representation on W 1 . 

For p G P the orthogonal projector I — p is denoted by p L . Let T s {p3f) denotes the symmetric Fock 
space over pS£ '. The separation of variables in finite dimensions extends to general symmetric Fock 
spaces owing to the canonical isomorphism of Fock spaces 

T p : JT = T s {2?) -> r s (p3T) ®T s {p x Sf), (23) 

for any finite dimensional projector peP, with T p Q. = QP^ ®QP ^ when QP^ and Q. p±3f are the 
vacuum vectors of the corresponding Fock spaces. We will omit the notation T p and identify directly 
the tensor products. 

Fix p G P. The tensor decomposition of the Weyl quantization comes from the Weyl relation which 
implies 

w^ + ^)=w^)w^')=w p ^)0W ] A^) 

for any G p3f x p L 3£ '. The symbols W p stands for the Weyl operator defined on the Fock space 

T s (p3f) and the Weyl quantization of b G ^(F), for any finite dimensional complex subspace F of Sf, 
is denoted by b^ eyI . Hence the Weyl quantization of b G y cy i(Sf) based on pSf equals 

b Weyl = I m[z)W{ V27lz) L p {d Z ) = ®^*) ■ 

In order to apply directly the finite dimensional results on Weyl quantization, we need to specify the 
correspondence of representations. 
On R d the Weyl quantization is often introduced as 

. {x-y).f .x+y , . d^dy 



b w <*{x,KD x )u{x) = f e^bi^^uly) 



(2nhy 



By a simple conjugation with a dilatation, it becomes a Weyl (Vhx,y/hD x ) where the position (x) and 
frequency (£) variables play the same role. An equivalent definition can be given with the help of the 
phase translations : 

T (loA) = e «^^ = (^ a »«) ' , K^u]{x)=S^I\(x-xo). 

It reads 

b Weyl {Vhx,VhD x ) = f ^\b\{y^)e lin ^ hx)+ ^ hD ^ dydr] 

JT*R d 

The symplectic form [[ , ]] and the scalar product ( , ) on T*R d are defined according to 

IfoO.Cfc*?)] =Z-y-x-n = -lm{x + i£,y + ir}) = -o(x + i^,y + iri) 
((x,^),(y,V)) =x.y + %.ri =Re(x + i%,y + iri) =S(x + i%,y + iri). 



After noting 



\[hx + Vhd x , Vhx — Vhd x 



2h. 
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the correspondence with the definition (1211 ) is summarized in the next table 

p% ~ C d T*R d 

(zuz 2 ) =S(zi,z 2 )+io(zuZ 2 ) z = e i6 (x + it,) {{x u E, y ) , (x 2 ,E, 2 )) = E, i .E l2 +x l .x 2 = S{z\,z 2 ) 

[[(xi,£l),(x 2 ,£>)]] =£,l-X2-X\.E,2 = -0{ZUZ 2 ) 
d d 

a{z)=a(Y J CCjej) a(z) = £ Wj(Vhd Xj + \/fac ; ) 

7=1 7=1 
d d 

a*(z)=a*(^ a^j ) a* (z) = £ dj (- V/j^. + Vfe; ) 

7=1 7=1 

[a(zi),a*(z 2 )] = e(zi,z 2 ) £ = 2/j [a(zi),a*(z 2 )] = 2ft(zi ,z 2 ) 

<£>(zo) = —7=(a(z )+a*(zo)) z = *o + ^o V2h(x .x + % .D X ) 
v2 

£( Zo )=lV(^ Zo )Q ¥ Co /.V () T k) (K- d '\-^) 
£o"> ko| = 1 Hermite function 

&eN fceN 

Once this is fixed, the general results on the semiclassical Weyl-Hormander pseudodifferential cal- 
culus ( IBoLell BoChl IIHeNill llHorl llMarl IINaNill llRobl ) can be applied for any fixed p e P. The notion 
of slow and temperate metric and weight depend only on the symplectic structure which is given by 
g(zi,z 2 ) = Im(zi , z 2 ). With such a metric the gain function X is given on p3f by 

xHz)- « with S?(T)= sup m! = sup 



rep;r\{0} g z (r) 5e P ^\{o} Si 5 ) 56 P ^\{o} 

With a slow and temperate metric g and a slow and temperate weight m, is associated a symbol class 
usually denoted S(m,g). 

After writing X = (x, € r*M d for the complete phase-space variable, the differential operator Dx is 
(D X ,D^ ) = (i -1 ok, i -1 ^ ). In the composition formula of symbols, the differential operator ^ [[Afp-D;^]] 
appears. After recalling 

1 ,„ .„ s . 1 



it equals 



^=2^ + ^) and d z = -(V x -N^] 



ih E 

-lD Xl ,D X2 ]] = -(d zr d- Z2 -d ?d .d z , 



We refer to [NaNi] for an explicit semiclassical writing of the Weyl-Hormander calculus within the 
Bony-Lerner presentation ([BoLe]) and with a general version of the Beals criterion following Bony- 
Chemin fllBoChll ) . 



Proposition 3.2 Let g be a slow and temperate metric on p^, dime (p3f) = d and let m\ and m 2 be two 

Neyl 



slow and temperate weights for g. For be 6 S p ^(m£,g),£ =1,2, the operator bJ e J^. acts continuously 



17 



on n keN D((N p s) k ) and on U keN D((N p ^ k ■ 
The symbol bi# e ^ 2 b 2 °fb\'p%' 



Weyl ,Weyl .. c 

o b 2 * s satisfies 



M d ^2-^^2)b l ( Zl )b 2 (z 2 ) 



l re 



\=Z2=Z 

j 



I -rr {^(d zr d- Z2 -d- zr d Z2 )) 6i(zi)*2(z 2 ; 

0<j<vJ- KZ ' 



+ s v R v (bi,b 2 ;e) 



where R v (bi,b 2 ;e) is uniformly bounded w.r.t e in the Frechet space S p ^ ( , j 
Vaillancourt theorem 



The Calderon- 



and the Gdrding inequality 



■ Weyl 
7 pJ° 



(&>0) 



> 



-C'p' kd {b)e 



respectively for b G <S p j»(l,g) ft G S p ^(X,g) . The index kdfor the seminorms p kd and p' k ) recalls 
the dimension dependent number of derivatives required in the estimates. 

The typical example Hormander metrics, which will be used here, are \dz\ 2 = dx 2 + d^ 2 (A(z) = 1) 

1 + \z\ 2 ) ■ Both of them split up in the (x,£) coordinates and the 



and 



\dz\' 

w 



dx 2 



;,c)> 2 + ((5)> 2 



Beals criterion of Bony-Chemin IIBoChl translated in the semiclassical case in [NaNi]- Appendix- A can 
be applied. Following the method recalled in [HeNi]-Chapter-4, this allows to check that functions of 
fully elliptic self-adjoint pseudodifferential operators are pseudodifferential operators, with an explicit 
knowledge of their principal symbol. In particular, this can be applied with 1 + +N p! t = (1 + 

|z| 2 )^J^ while noticing that 1 + ed ™ p +N p ^ is a fully elliptic operator in S((z.) 2 , ^jr) (such a result 

with £ = 1 can be found also in [Hell]]). 



Proposition 3.3 Fix p G P, fix the exponent j£l and let N p ^ = dT(I p g) be the number operator on 
r s (p2f). For any s G R, (1 + ^ + N p s) s/2 can be written (b(s,e))™%! with e- l (b(z;s,e) - (z) s ) 

uniformly bounded in S((z) s ~ 2 , yrr) • 



3.3 Weyl quantization and Laguerre connection 

In this paragraph, the relationship between the Wick and Weyl calculus is checked in the infinite dimen- 
sional setting. It specifies the relation between the representation of the Weyl algebra, generated by the 
W(%), and the number representation which puts the stress on Wick symbols or Hermite states z . This 
relies on the introduction of Hermite and Laguerre polynomials, recalled below. 
Let h n [x) denote, for any n G N, the n-th Hermite polynomial in C: 

jn [i/2] | 

Those classical polynomials are also given by the generating function 
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Lemma 3.4 

(i) For any the following identity holds in r J%?fi n : 



,- / r- \ Wick 

1^1" , (iy/2S(Z,z)\ 



2»n\ n \ |^| 
( ii) For any n, j,k £ N the estimate 



l m (N)oh n (iV2S^,z)) ol {ke} (N) 



Wick 



n\ 



<(l + 2 v /2(* + j)e|«|)" _ 



holds for any E, £ 3?. 



Proof. Using the generating function (T23T ) with t = and x = implies the equality of the 

Wick symbols 

to 2 "" ! "V l^| )' 

Nevertheless the equality of the the series of Wick quantized operators has to be checked. 
Recall that elements of J#fi„ are analytic vectors with infinite radius of convergence for the field opera- 
tors. Hence the sum 

oo -ft 

n=0 ' 

is absolutely convergent for all € 3f. Therefore to prove (£) it is enough to compute the Wick symbol 
of <I>(^ )" for all n. Indeed using the Wick ordering rules, we have 

[n/2] , I e |2r n-2r 

~ V2V(n-2r)! 2 r U " 2r VS7 VS; 



= ^ £r!(n-2r)! ^ ^WP 7 £ C "- 2 ^' ^ ' z) " 

J^T ^ n\ ( 2V2S(S,Z) ,n-2r\ 

2" \hr\{n-2r)\ K |§| ' J ' 

To prove the second statement (ii), take y/jt £ V* an d V; ^ V 7 ^° an d write 

Wit* n ! / \n-2^ W,dt 



(V0-,^(i^5(^z)) y t ) = £ (n _ 2r) , r , (VO, ((2iV2^,z))"' 2r ) i^). 
Using Lemma [2~4l one obtains 

(^,/j, I ^V25( l §,z)J y/*) < |VOi V J^ly*lv*^L ( w _2r)!r! 2 (^ + j)^ 1^1 

^ vol 

(l+2v/2(^ + 7)£|§|r r — - 



in— 2r 



[,/2]! 
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□ 



The Laguerre polynomials are defined by the formula 



4°(0 = E(- 1 )" 



(k+j)\ 



f\ t £ C. 



m=0 {k-m)\(j + m)\m\ 

The following proposition gives the Laguerre connection (see IFoll .fRipO. 



Proposition 3.5 For z,% £ ^ with \z\ = 1, the next equalities hold according to the ordering of j and 



nyle 



(iy-j^Lf- j \\(l;,z)m,z) k - j e-\S\ 2 / 2 i£k>j, 
[ {^- k ^t k \\^,z)\ 1 ){z^) j - k e-^^ i£j>k. 



(26) 



Proof. Let us establish the expression of r f[z® k ,z®^] in the case k > j. The case j < k is similar. Using 
Lemma [33 one obtains 



y[z»*z«y](-*=) = (z®i,WU-J)z® k ) 
K\/2e V e 



V 



Now let use the explicit form of h n and Proposition [23] We obtain for |z| = 1, 

[n/2] n -2r i"\^\ 2r 



Wick 

y [z ® k ,z® j ](^=) = E E E *y ^ c° n _ 2r e r -i (z^,{(^.yur 2r - s ) z^ 



^V2e to to to 2 r r!(n-2r)! 



[n/2]»-2r ;n|,£|2r 



E E E or, 



j n l£l 



( — l)*Jfc! 



V-toto 2'r! s\(k-j + s)\(j-s)\ 
The last term gives the claimed identity. 



□ 



3.4 Anti-Wick Operators 

The Anti-Wick quantization is introduced by a separation of variables process like the Weyl quantiza- 
tion. For a given p £ P, set p = 1 — p, and use the tensor decomposition ([23l . The Weyl operators 
on p% and p^2? are denoted by W p (£i) and W p ±(&) with ® L = Wp(^i) ® W p ± (&) . For any 

£ £ pJT, the coherent state E p (£,) is defined by E p (t;) = W P (^-)QP^ . Introduce the projector on 
after tensorization with Ir s tp± 

piT 9 § -if = (\E p (t;))(E p m®Ir s ( P ^) ■ 
The Anti-Wick operator associated with a symbol b € ^" c> ./(iF) based on is then defined by 



lA-Wick 
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The above formula can be first considered in a weak sense or as a Bochner integral when i 6 y(pf ) 
and the bounded projector P| is continuous w.r.t. The finite dimensional identification of the Weyl 

symbol of \W p (^-)Q. p ^)(W p (^-)Q. pS |, can be deduced after completing the table of correspon- 
dences in Subsection P 



p2f~C d z = x + i% T*R d 

r s (p3?) ~ r,(C d ) , e = 2h L 2 {R d ) 

Ma) =w P C§zo)ar* ? = u 

\a pS ") (£l p3£ \ = y Weyl {K)- dl2 e-'i- y i = g Weyl (Vhx, VhD x 

y{z)=2 d e 4= withg(x,£)=2 d e — ~ 
From the conjugation 

^ ^a Weyl (Vhx,VhD x )T* ^ =a(.-xo,.-So) Weyl (Vhx,VhD x ) 

the above correspondence gives 



E p ^))(E p ^)\ = Y f eyl with J^z)=2 d e 



— 



Hence the usual finite dimensional relation between the Weyl and Anti-Wick quantization now reads 
(after tensorization with I^fpi %\) 



/ r.i2 \ Weyl 



I -A -Wick 



(27) 



P s (ne /2) &mpS 

\ 

f &[b]£) W(V27l^ e-^'LpidZ), (28) 



for any b G ^(p^f) by setting 



>*Y(z)= b(z)r(z-z')L p (dz l ). 

p3f J p 2? 



From (1271) . the Anti-Wick quantization can be extended to symbols in 5(1, \dz\ 2 ) with the next properties 
(see IfHMRl). 

Proposition 3.6 Fix p 6 P. Let b G S p ^(l, \dz\ ), the following statements hold true: 

(i) Ifb > then b A - Wick > 0. 

(ii) i* 4 -***!^ < 

C »z) 77ie comparison with the Weyl quantization is given by (1271 ) with the estimate 

\ b A-Wick_ b Weyl^^ £CdpMe 

where the constant Cd > and the seminorm pk d depend essentially on the dimension d = dimpSF. 

A variation of it holds when b € (^b(p3f)), when ^b(p^) denotes the set of bounded (Radon) 
measures on p3f and comes directly from (|28l ). 



21 



Proposition 3.7 For any pGP and any b G 1 

asymptotically the same: 

X{m \ b A-Wick_ b Weyl 



(p3?)), the Anti-Wick and Weyl observables are 
0. 



^o 1 

Proof. Recall that b We * 1 can be defined for any b G S"'(p3f) as a continuous operator from r\ke¥\D{N k p ^) 
~ ^(P^) to U teN D(A^ J »)* ~ ^'(P^), with d = dim/^ and (EH) is still valid for such a symbol. 
Assume ^7? = V G J(b{p^)- The identity 

(V, (& Wey ' - ft* - ***)?) = ^ W(V2^)<p) ^l-e-^l«l 2 ^ rfv(§) 

holds for any <p, y/ G rikeND(N k ^-). This implies 

_yb.-Wick\ 



< 



□ 



3.5 Weyl quantization and specific Wick observables 

In finite dimension, that is for any fixed p G P, polynomially bounded symbols can be introduced af- 

I , |2 

ter considering the class of symbols U se M.S p s{(z)'\gp) where g p is either the metric \dz\ 2 or vpr on 

According to Proposition 13.21 it is an algebra with the Moyal product, # e / 2 , associated with the 
composition of Weyl quantized observable with a complete asymptotic expansion of oi# e / 2 Z?2- For any 
m,q G N, the finite dimensional space & m ^(jj3?) of (m,g)-homogeneous polynomials on 3? is con- 
tained in yps((z} m+q \Sp)- The comparison between the Weyl and Wick quantizations is symmetric to 
(ETJ) (seeEeSh)): 

Wick 



( 



V»g 



Weyl 
p3f 



E/2 



\ 



P 3f (7T£/2) dim ^ 



For polynomials the deconvolution is possible and we get for any m,q G N and any b G 



where the symbol c(e) equals 



c £ 



/: 

V 



1 pi 



* 



(we/2) dim ^j 



and is uniformly bounded in S p % -{{z) m+q ~ 2 ,g p ) w.r.t £ G (0,e). 

The space ^ m ^ q {p^) is identified with a subspace of & mA {3?) and even of any g? r m q {!¥) for any 
r G [1, +°°] with 

^be^> m , q {p^), VzG 3T, &(*) =fe(pz + p ± z) =fe(pz) 

5 = / ,®«o5o P ®'« = r,(p)5r,(p). 

After tensoring the finite dimensional comparison with /^(p 1 iT)> we nave P rove d 

Proposition 3.8 For any p G P, a«j m,o G N, f/ie class of symbols ^ mc/ (p^) is contained in q {2f) 
^Sp^{{z) m+q i§p)- Moreover the operator e~ 1 (b Wlck —b Weyl ) can be written c^^ 1 with c E uniformly 
bounded in S p %-((z} m+9 ~ 2 ,gp) w.r.t e G (0,£~). ( The metric g p can be either \dz\ 2 or $L on p2f.) 
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4 Coherent and product states 

We distinguish the coherent states E(z) = W(^z)Q (resp. the projector \E(z))(E(z)\) from the prod- 
uct or Hermite state z® k (resp. the projector \z® k ){z® k \). Although they are intimately related, the 
asymptotics of coherent state E(z) tested on Wick, Weyl or Anti-Wick observables encoded exactly the 
geometry of the phase-space 3?, while the asymptotics of the product states z® k , ke —> keeps track of 
the gauge invariance 

V0g[O,2tt], \{e w zf k )((e ie zf k \ = \z® k )(z® k \ 
with variations according to the quantization. 

Proposition 4.1 Fix z, G 3£ with \z\ = 1. 
(i) The convergence 



lim r[z®*,z®*- M ](§) = ^- r e 2niS ^^e- ime dB 
e^o 2n Jo 



holds for any fixed m G N by setting z 6 = e lS z. 
(ii) The coherent state E(z) = W($z)& satisfies 



r[E(z),E(z)]^)=e 2niS ^e-^ ^V*^ 

5' 



Proof. /) Set j = k — m and compute ^[z® jZ®- 7 ]^) with ^ = -JU according to Proposition 

y [z « iZ « K _^ ) = (ir yZ t M(£ K5 ' jZ) |2 )( £)»/2 (? ', z) - e - e |5'P/4 

The bounds (e£) < C and £~ =0 < °° imply 

lim = (o* £ ,J~ 1)J ■ l Kr,a)i 2, (r,z) w , 

V27T 32T"j!j+ffl| 



fe -> 1 



owing to Lebesgue's theorem. A simple series expansion e' = £~ =0 ^ for ? = iV2S(z. e , £') gives 

i- / ^(^V^dO = (f )» £ l(r,z)l 2 -(r,z) m - 
27T Jo S2"J"j|(i+»l)! 

h') is a straightforward consequence of (f20l >. □ 
The next result specifies the similarity and the differences between the product states and the coherent 
states in the mean-field or semiclassical limit. 

Theorem 4.2 Let z G 2f and m G N be fixed with \z\ = 1 arcc? set z d = e l0 zfor 6 G [0,2jt]. The next 
limits exist as e — > 0, ke — > 1. 
W For ft €^(2"), 

1 r 2 ' 1 

lim (z»*- m ,fc w ^z»*> = lim l z ®k-m^-Wick z ®k\ = _ b(z 9 )e- im9 de. 

e-,0 e-,0 2?T , 

*e — 1 *e -. 1 
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Meanwhile the coherent state E(z) satisfies 

lim(£(z) , b Weyl E(z)) = lim(£(z) , b A - Wick E(z)) =b(z). 

e^O e^O 
(ii) For be 3P p>q (3f), with p,qe N fixed, 

1 f^K 

lim (z®*- w ) & w *z®*) = 5„_„ m = / b^e^dd . 
e^o r 2n Jo 

*e — > 1 

Meanwhile the coherent state E{z) satisfies 

V£>0, (£(z),ft w '' c *£(z)>=ft(z). 

Proof. Set j = k — m, with m G N fixed. 

For (i), fix ft G ^/(JQ. The definition of 6^ in 42}, says 

(pj^weyi^ = f ^lb]^)( z ®i, W {V2nt;)z® k )L p {dt;) 

JpS 



ips 



Since J^"[ft] G ^(p3f) and ^[z®*,^]^) converges pointwise according to Proposition ^. 1 1 Lebesgue's 
theorem yields 

lim ( z ®j,b Weyl z m ) = I &\b]{£)(^- r e a%s ^^e- ime de] L p (d$) 
«-»<> JpS' \2n Jo J 

ke -> 1 

1 r 2n 

= si *V-<«- 

The limit with Anti-Wick observables is a consequence of the formula (|28T ): 

Jpsr 

The statement about the coherent state E(z) is even simpler by referring to Proposition 14. 1 1 (ii) . 
Let us prove (ii). The statement (ii) of Proposition [23] gives 

/ ®j iMick ®k\ _ c+ / 4 s ( 7 ®q U 7 ®P\ 



We conclude again with . / \ItSf 1 as * ^ °°- D 



5 An example of a dynamical mean-field limit 

In order to illustrate the general presentation, the standard example of the mean field derivation of the 
Hartree equation from the non relativistic Hamiltonian of bosons with a quartic interaction is considered. 
Two standard methods are considered: The coherent state method (see [Hep ][GiVel or MCasl for a rapid 
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presentation) also known as Hepp method and the propagation of chaos approach with a truncated Dyson 
expansion according lFg%llFKPl[EgYflESY2l|gpo|. 



Consider 2? = L^(R d ,dx) and take V G L^(R d ,dx) such that V(-x) = V(x). The polynomial 
Q(z) = ( z ® 2 , Q Z ® 2 ) is associated with the operator Q G if((g) 2 2f ) defined by 



m(xi)w(x 2 ) i-> -X2)m(xi)w(x 2 ). 

The Hamiltonian is defined as 

H e = dT(-A) + 0™* 

where —A is the Laplacian of W 1 , while denotes the free Hamiltonian dr(— A). It is well known 
that H e is a self-adjoint operator on Jt? (see MGiVell ) and the quantum time-evolution group is denoted 
by U e (t) = e~'e He while U°(t) = e~^ H ° = F(e itA ) stands for the free dynamics. Although the Wick 
quantization of non continuous polynomials has not been introduced here, this Hamiltonian appears as 
the Wick quantization of the energy functional 

h(z) = I |V Z | 2 dx + Q{z). (29) 
It is also well known that the mean field limit is the nonlinear dynamics issued from the Hartree equation 

id t zt = -Az t + V* \z t \ 2 z t = dzh{z t ) (30) 
with initial condition z,o = z 6 3? . 

An important property of the dynamical groups U E (t) and U®(t) is that they preserve the number 
U e (t)*NU e (t) =N, [H e ,N] = [H° 8 ,N] = [Q Wick ,N] =0. 



Remark 5.1 All the results of this section can be easily adapted with a self-adjoint operator A on 2f 
and a polynomial Q{z) G ©^^^^(iF). Nevertheless it is better to stick to this concrete presentation 
which fits better with a widely studied problem. 

5.1 Propagation of squeezed coherent states (Hepp method) 

In finite dimension it is nothing but checking the propagation of gaussian wave packets. Although it 
works only for some specific states it is a direct and very flexible method. Moreover it agrees very well 
with the phase-space geometric intuition. Extensions with more singular potentials or about the long 
time behaviour have been carried out in [Hep)IGi"Vel . 

Proposition 5.2 For any zo G 3f, the estimate 
holds with 

id t z.t = -Az.t + (V* \z t \ 2 )zt , z t =o = z.o (31) 

co(t)= f Q{z s )ds (32) 
Jo 

ied t U 2 (t,0) = [dT(-A) + Q 2 (t) Wkk p 2 (t,0) , £/ 2 (0,0) =/, (33) 

Q2(t,z) = \ MQ(z t ) ,z® 2 ) + (z 82 , diQ(z t )) + 2(z, d- z d z Q(z t )z)) , (34) 

(d 2 Q(z t ),z® 2 ) =2(Qz? 2 ,z® 2 ) G 

(z , dgd z Q(z,)z) =4(zVz t ,QzVz t )e & 1A (2f) . 



< Ce C|V|^{ Z0 ) 2 (|f| + l) gl/2 
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Proof. This proposition says that the evolution of a coherent state is well described after applying a time 
dependent (real) affine Bogoliubov transformation like the ones considered in Proposition 12.111 
It is sufficient that 



e i ^e i ^W(^z t )U 2 (t,0)a = e i ^r(e itA )e i ^W 

ie 

J~H E rf„itA\ 



V2 
ie 



remains close enough to Q.. The quantities U E (0,t) = e ! 't^r(e SA ), U 2 (t,0) = r(e- itA )U 2 {t,0) and 
z t = e~" A Zt solve the differential equations 



after setting 



ied t U e (0,t) = -U e (0,t)r(e- itA )Q Wick r(e itA ) = -U e (t,0)Q(t) Wick , 
ied t U 2 (t,0)=r(e- itA )Q 2 (t) Wkk r(e itA )U 2 (t,0) = Q 2 (t) Wick U 2 (t,0). 

idrz t =e- UA (V*\e itA z t \% lt % = d- z Q(t,tt) , z =z , 



Q(t,z) = Q(e ,tA z) and Q 2 (t,z) = Q 2 (t,e' tA z) 



(35) 
(36) 
(37) 

(38) 



The main properties of U 2 (t,0) are derived in MGiVel Proposition 4.1] and in particular we know that 
U 2 (t,0)Cl belongs to the domain of the closure of any b Wick with b £ ®t% e ^ P , q {^)- 
The differentiation of the Weyl relation Q on J£/ !H says 



V2 

ied t W{^Azt) 

l£ 



-Re{z t ,id t zt) + V2^(id t z t )]w( 



V2\ 



l£ 



-Zt) 



[-Re^,^^)>+a*(ay^(^))+fl(^ f (^))]W(^) 



-Re(z f , d-MM+Mz, dzQt{lt)) Wick ] Wi^zt) ■ 



The translation property (iii) of Proposition |2.9| then leads to 



e ^Hs e i^ W (2_z Zt )u 2 (t,0)Q. - Q. = - I U e (0,s)e l ^W(—z s )^(s) wlck U 2 (s,0)£l ds 
after testing both sides on J#fi n and setting 

*S(s,z) = -Q(s,z + z s ) - co'(s)+Re(z s ,d- z Q(s,z s ))+Re(z, ^<2,(z,)) +Qi(s,z) 
= -Q(s,z + z s ) + Q(z s ) + (z,d,QSs)) + (d z Q s (z s ) , z) + Q 2 (s,z) . 

The last equality is given by our choice of co(t) in (l32l) . It suffices to find a uniform estimate w.r.t 
s G [0, t] of the squared norm 



.V2 



1 



\Wickfr 



Wick,* 



\Wickfi 



e-W(s) mck u 2 (s,o)a\^ = fr 2 (n, u 2 (o,s)^(s) wlt% *^(s) mck u 2 (s,o)a) . 



(39) 



The important point is that the symbol £f (s) vanishes at the second order at z = 0. More precisely it can 
be written 



with t* M (s) e &> p , q (3f) 



and 



W P ,qi S ) 



^ <C PA \V\ Ir \zo 



\4~P-<i 
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Owing to Proposition 12.61 and Lemma 1231 the operator s& (s) Wtck >* s/ Wlck {s) takes the form 

2 

k=0 6-2k<p+q<% 
With \^k,p,q{ s )\g>(yi> eg \jq ag\ — Ck,p,q (zo) ■ 

The estimate of every term 

e k ~ 2 (n,t) 2 (0,s)^ k , p , q (s) Wick Ms,0)&) , P + q>6-2k 
is given by the Lemma 1531 below and yields the result. □ 



Lemma 5.3 Consider the time dependent Wick operator Q 2 defined by d34|) d38b and parametrized by 
Zo G 3f. Consider the associated unitary operator U 2 (s, 0) defined by (136b - For any p, q € N, there exists 



a constant C pq such that the estimate 



p+q 

ly(V"f,V"5") £ 2 



/ Q,U 2 (0,s)b Wick U 2 (s,0)O) \ < C M e c P^M 2 (\^) \~ b \ 
holds for any b G £^ p ^(3f) and any s£K, 

Proof. By introducing an anti-unitary operator Jz. = z. The M-linear operator o\Q 2 (t) can be written 

o\Q 2 (t)z = R{t)z + R 2 {t)z. 

The definitions (1341x1381 ) ensure that R{t) is a bounded operator strongly continuous with respect to t G 1L 
and that R 2 {t) is a Hilbert-Schmidt operator which depends continuously on t G E in the Hilbert-Schmidt 
norm. Moreover the following uniform estimates hold 

l^(0l^)<2|V| i .|z | 2 , \Ri{t)\^ m <2\V\ L „\zo\ 2 ■ 

Hence the equation 

id t $ 2 = b\Q 2 (t)<t> 2 =R(t)<t> 2 +R 2 (t)J<t> 2 
defines a dynamical system of bounded M-linear operators with the estimate 



More precisely the Duhamel formula 

r h - 



3> 2 {t 2 ,ti) = Te 



i j 1 Te- i5 ' lR{s) ds R 2 (t)M 2 (t,h) dt 
Jti 



implies that the R-linear operator & 2 (t 2 ,h) can be written 

<& 2 {t 2 ,t { ) = B(t 2 ,t\) +B 2 {t 2l h)J 
with \B(t 2 ,h)\ n2r) + \B 2 (t2,h)\je2 (T) <C\V\ L „ \z \ 2 (\h-h\ + lJgCfe-ill^-NI 2 . 

According to Proposition 12. 1 1 1 for any c G ®p+g= m ^p,g(^°) and any t G K, the polynomial c(t,z) 
c(<J> 2 (0,0z) belongs to @ p + q = m @ > P ,q(2?) with 

p+q=m p+q=m 
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Applying the characteristic method, that is differentiating c(z) = c(t,<t> 2 (t,0)z), shows that c(z,t) solves 
the equation 

id t c{t,z)+d z c(t,z).d I Q 2 {t,z)-d z Q 2 (t,z)d z c{t,z) =0. 

Thanks to the Wick calculus in Proposition |2T6l and the fact that U 2 (t, 0)Q. G Pi ke n@(N k ) (see MGiVel 
Proposition 4.1]), this leads to 

id t U 2 (0,t)c(t) Wick U 2 (t,0)D. = U 2 (0,t) (e-\c Wick (t),Q 2 (t) Wick ]+id t c(t) Wick ) U 2 (t,0)Q. 



= 2 (o,O§ ({c(t),Q 2 (t)} {2) ) W ' ck u 2 (t,o)a- 

Take b G ®p+q= m &p,q{3?) an d apply this result with c defined by c(s,z) = b(z), which means 

c(O 2 (0,j)z) = c(s,z) = b(z) 
or c(z) = fe(*2(j,0)z) G ®p^=mo^(30 



With £ I^WI^^T)^^ 1 ^ 1 ^' 1 '^ I |W)| W , lV ,fl 



'1 <4,Ml~W 2 (|*|+1) 
'mo 

p+q=m p+q=m 

This leads to 



a, ^2(0,^^2(^0)0^ = ^a J c^a)+y'^n,a f ^2(o,Oc(0 W! ' c ^2(f,o))aJ) & 



By noticing that the symbol {c(t),Q 2 (t)} vanishes when mo < 2 or belongs to ® p + q = m -2&p,q(3f) 
with 



£ a|^{c(o,^(/)} (2) | <c|y^|z | 2 I |W)| W5iVIJ) 



if ( V p ^,V ? 



p+q=mo 

the result is proved by induction on m and by using x" < «!e* for x > 0. □ 
5.2 Truncated Dyson expansion 

We focus now on the propagation of chaos point of view which has been considered by several authors 
in IESY1I |ESY2) II B GGMll MFGS I . In the bosonic setting Hermite states tested on some Wick observable 
is exactly the BBGKY hierarchy. For example the reduced one particle density matrix can be defined 
as TrfpiA] = Tr[pdT(A)] = Tr[p£? Wlck ] with s£ (z) = (z.,Az.) . While reproducing the Dyson expansion 
analysis of JFGS], we check here that a full asymptotic expansion can be written, when Wick observ- 
ables are tested after the suitable number truncation. 

The strategy of the proof in [FGS ] relies on an analysis of the Schwinger-Dyson expansion of a time 
evolved observable U e {t)*ff U e (t) is given by 

U e (ty0U e (t) = t + j Q dh--J o dt n [Ql ick ,--[Qf c \ff t \---\ (40) 

where ff t = U°(t)*0 Uf}(t), Qj kk = U°(s)*Q Wick Ujfrs). The commutation relation in Proposition 
(iii) yields 

QY ick = ({(e hA zf 2 ,Q(e isA zr 2 )) mck , 
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or shortly Q s (z) = Q{e' z) and we shall set more generally for b G & p . q (3f) and s G 



ft, G : Vz G 3f, b s {z) = b{e ls \) . 

Although the convergence of the series can be proved as an operator acting on V 2f, with k G N fixed, 
the e-asymptotic analysis is done with its truncated version 

U e (t)* G U e {t) = ff t + f^{ i -r fdtr- f n ~ 1 dt n [QZ ic \---[QY 1 ick ,0 t }---) 

£ x e Jo Jo 



, -y I dti- I d ti Ue(t e yu2(t t )\QZ K *,•••[&?' "^■■•mtyueiti). (4i) 

£ Jo Jo 

The Poisson brackets analogue of the multicommutators will be necessary. 

Definition 5.4 For n, r G N, r < n and any fixed b G & PA (3f), the polynomial ci"^ (t\ , . . . ,t n ) is defined 
by 

&\t„,... A ,t) = l- r £ {Q tn ,---,{Q tl ,b t } {ei) ---} {£n} ^^p-r+n, q -r+n{2f), (42) 



e/6{l,2} 



and Cr(t\,...,t n ,t,z) denotes its values at z G 2f while c[ n \t\, . . . ,t n ,t) or simply Cr denotes the 
associated operator according to Definition 12. 1 1 . 

We shall prove. 

Theorem 5.5 Fix p,q G N and assume b G & Ptq {3f). Then the asymptotic expansion 

Wick 



u e (tyb Wick u e (t) = J £ £ rJ £i n dt 



r=0 n=0 







( dt\-- 


■ / dt n 


Jo 


Jo 



Cr"\t n ,-- ■ ,t\,t) 



+ £ l Rl{£,t) 



holds for any £ G N and any 8 > in (\J k 3f, \J k ~ p+q j§?) with the uniform estimate 

\Rl(£it)\^ryk ^yk-p+q ^ < when ke < 1 + 5/2 and A{\ + 28)\t\\V\ Lr <\ . 

A particular case takes a more explicit form. 

Theorem 5.6 Take b G B^ p ^ q {2f). Let z G 2f be such that |z| = 1 and call z t the solution to A30i with 
Zo =z. 

( i) Then the expansion 



{z® k - m ,U e (tyb Wick U e (t)z 0k ) = 8, 



p—q,m 



£V p("\t,z,k,e) + Ot(e e ) 

r=0 



(43) 



holds as e — > 0, ke — > 1 by setting 
P m (t,z,k,e) = Hz,), 

<*%«,*,«> = t lV N/t!(t " m,!£ " +,+ * 



(k- (p + n-r))\ Jo 



t n -\ 



dh--- I dt n ci n \t n ,--- ,t h t;z), (44) 



and as soon as 4\t\ < 1 . 
( ii) More generally the limit 

lim (z® k - m ,U e (tyb Wick U E (t)z® k ) = 8 p ^ m b(z t ) 



ke -> I 



holds for all times t G 
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Corollary 5.7 In the specific case m = 0, q = p, the expansion (1431) takes the form 

(z® k , U e (t)*b Wick U e (t)z m ) 



£e*£f / dt x - / d*„ 



[£of^(te)C^(r B ,.--,ri,j; Z )] + 0(e') ) 

7=0 

where the coefficients OCj" ( K") are polynomials in K given by 

p+n—i — 1 

£ a^"(K)e J = k(k- e)(K-2e) ■ ■ ■ (k- (p+n-r- l)e), 

;=o 

a?i<i convention that otj" = w/ie« j > (p + n — r) or r > n. 

Proof. We are considering the particular case p = q, m = 0. Setting k = ke = (k — m)e gives: 

(k-\p + n-r))\ =K{K - £){K - 2£) - {K - {P + n - r - m - 

Putting together the terms of order e s , s less than I — 1 in Thm. 15.5t h). yields the result. □ 
Before proving Theorem [53] and Theorem [531 l et us collect some technical preliminaries. 

Lemma 5.8 Forb G ^ p , q { 3f) the identity 

1 n / , n \ Wick 

y n [Q^ c \-\Ql ic \bT ck \\ = 1^ (&\t n ,...,t h t)) , 

/ioWs w/f/i ?/ze symbols C^fi {t\ , • • • defined according to (|42l) m Definition 15.41 
Proof. Proposition |2.6l provides the induction formula 

^^{e^^^ + ^G^C^^^, (45) 

with c£ = if / < r or r < 0. In particular, we get 

C { n) = {Q,„,-, {QtM}. 
A simple iteration of (|45| > yields the result. □ 

Lemma 5.9 Let b belong to ^ l Piq (3f). 
( i) The estimate 

Sl ^(v p+I ir,v 9+1 ir) - + ^ ' V ' L °° ^(V'^.V**)' 

toto fry M ftmg Si = {p l }) , {q l x) , d! +l d! +1 {QsM W G ^(V p+1 V* +1 
( h'J Similarly, the inequality 



-2 



( Hi) For any n£N arcc? r G {0, 1 , . . . , n}, the operator Cr associated with the symbol (t n ,...,t\,t) G 
& 1 p+n-r,q+n-r{^) according to Definition \5.4\ satisfies 



c {n) 



<2 C n {p + n-r) — — — — - — \y\L"V>W(\i'ar,\nar)- 

/p+n-r oy> \/q+n—r oy\ \j) 1 ) . 



when p >q with a similar expression when q > p ( replace (p + n — r,p—\) with (q + n — r,q—\)). 
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Proof. The statements (i) and (ii) are particular cases of Lemma 12.51 The estimate in (iii) is a conse- 
quence of (i)(ii) and the definition (l42l) . □ 
Proof of Theorem |53J Set j = k - p + q. Since U e (t) and U°(t) preserve the number like the 
equality 



U e (t)*b wu *U e (t) 



n=0 



tn-1 



dt,[Q7j a ,--[QZ'"',bY>*\---\ 



dt t u e (t e yu°(t e )[Q™ ck , ■ ■ ■ [Ql ic \bT ck ] ■ -]u°(t e yu e (t e ) . 



derived from gB holds in =Sf (V* 2", V 7 Then Lemma EH implies 



-l 



u E {t)*b Wick u e {s) = £/" / rfr 



n ■ ■ • 



t„-l 







r / <i? 



i • • • 



n=0 1/0 1/0 r=0 

'f- 1 n « T M 1 Wick n 

dt t U e (t e )*U°(t t )e* Cf{h, ■ ■ • ,fi,f) U%(t t )*U e (t t ) 



(46) 
(47) 



£-1 

^C/ e (^)*C/ e °(^)£e r C^(fc,-,*i,0 ""~U°{t t )*U e {t t ). (48) 

r=0 



Keep untouched the part (I46i>-(|47T) and iterate the Dyson series on the third term (l48l) . While doing so, 
use the formula 



n Wick i-\ 



r=0 



Cr"\tn, ■■■ ,tl,t) 



Wick r -iWYefc 



r=0 



+ - 



{ft„ +1 ,cW(wi,---,?i,0} (2) 



(49) 



inductively for n = £ + 1 , • • ■ , M — 1 . After M — I steps, collecting the factors of e e yields 

M-l -f 

U e (t)*b Wki U e {t) = £ f / • 



£i n / dh 



-i M I dh 



M-l r f min(f-l,n) n ^ 

£ £ r a n \t n ,--- ,t u t) 



tld-\ 



dt M U e (t M )*U°(t M ) J £e r C ( r M) (t M ,--- ,t u t) ' "'~U°(t M yu e (t M ). (52) 

r=0 



Assume that for 5 > there exists a constant Cg such that 



(50) 



dt n u e {t n yu Q e (t n )- [{& n 4-i ('»-!>••• ,*i,0} (2) J u Q e (t n yu e (t n ) m 





ft,,- 1 


j dh-- 


■ / dt n 


Jo 


Jo 



Cr {tni ' ' ' )t\ it i 



<c s . 



(53) 



According to Lemma[231 the first term d50]) of d50l)(l5Tl)(l52l provides in U e (t)*b Wick U e (t) | v * s the partial 
sum of a convergent series in J^{\/ k ^,\J k ~ p+q when ke < 1 + | . With the same argument the 
remainder term (l52l) vanishes as M — > oo and < 1 + f . By referring to Lemma [5T91 (ii) and again to 
Lemma 1241 the factor of e in (f5TJ) is associated with a series which converges in 
as M — > oo uniformly w.r.t. (A:, e) when < 1 + f . The sum of the series is simply denoted by R((t,e). 
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Let us prove (l53l to finish the proof of (ii). Lemma l2~4l and Lemma [5791 say 



r=0 





rt„-i 


fdh-- 


• / dt n 


h 


Jo 



Cr ^ {t n j * * ' j t\ , t ) 



\t"\ 



<Y(l + 8) n Y l -t max 



r=0 



2 n-r\ t n 



i?(V" + "~' 2r,\/« +n -' '%) 

(p + n — r—\)\ 



< 1(1 + 5)"! ^C:A(p + n-r)(p + n-r-l)Y K1 ^ >- |V|£- |^ (V ^, V ^) 

r=0 V" 1 /" 



2 n- 



OO 

<2P^(\ + 8)H n \t\'\n + p) 2 '\V\l.\b\^,^ N ^ ) . 

n=i 

The last r.h.s. is finite whenever \\t\ \V\ L ~ < + The condition (1 + 28)4\t\\V\ L ~ < 1 is sufficient 
and provides the uniform bound Cg in (l53l . □ 
Proof of Theorem 15.6b Set j = k — m. By Theorem 15.51 the right-hand side of (l43l) vanishes when 
m 7^ p — q and the convergence of the series in 

JSf (V* 3T, \J k -P +c i 2?) combined with Proposition E/J-ii) 

implies 

(z® j ,U e (t)*b Wick U e (t)z® k ) 

l-\ co / 



y £ ry - / fe!j!e/^-2("-r) 

rtn-i 



[ dh- - [" 1 dtndPlfn,— ,h,t;z) + S (s e ), 

Jo Jo 



when ke < 1 + |, for any 5 > 0. By considering the limit e — > 0, he — > 1 every factor 



k\j\eP+i+ 2 (. n - r ) 



(k-{p + n-r))\(j-{q + n-r))\ 



converges to 1. Therefore this proves (ii) for small times t such that 4|f||v|z,°° < 1 up to the identification 
of the first term as b(zt)- From our definitions we know 



b( Zt ) = (z? q M? p )=b t (e-^z s )\ s=t . 

By setting w s = e~' sA z s , the quantity b{z t ) equals 

b{z t ) = b,(w Q ) + / d s [b t (w s )] ds = b t (w ) + / d s w s .d^b t (w s ) + d z b t (w s ).d s w s ds 
Jo Jo 

Moreover the equation d30l) has the equivalent form with the vector w s = e~' sA z, s and W s 
id s w s = e-^Qiz,) = d- z Q s (w s ) - id s W s = d z Q s (w s ) . 

Hence we get 



b(z t )=b(w ) + i [ {Q tl ,b,}(w tl ) dt y . 
Jo 
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An induction with wq = z and the convergence of the series already checked yields 

°° [' t Pt — 1 

KZt) = Z / d h— I dt n C^\t n ,...,h,f,z). 

~ Jo Jo 



Now let us prove the limit (i) for all times by following the argument in MFGSI1 , | |Spo| . Assume that 
the result is true for \t\ < ■ Let s be such that \s\ < \/4\V\l-»- The convergence of the series given 
in Theorem [53] and the fact that U e (t) preserves the number gives 

(z® j ,U e (t + s)*b Wick U e (t + s)z® k } 

OO ft pg pg _j 

= £r" £V / d Sl - / ~ ds n (z^,U e (t)*[d n \s n ,---,s 1 ,s)} Wick U e (t)z^) 

n=0 r=0 jQ jQ 

OO p$ _j 

= Yf / dsi- / ds n (z 0j ,U e (tr[C^(s n ,---,s h s)} Wick U e (t)z 0k ) + O s (e) (54) 

n=0 JO JO 

with an absolutely and uniformly convergent series in the d54l when ke is close to 1. Hence the limit 
e — ► 0, ek — > 1 and the sum £"=0 in (1541 can be interchanged when 4 1 ^ 1 1 V | < 1. An induction on 
K = 0,1,2... finishes the proof. □ 

5.3 Coherent states and Wick observables 

We show here that information on the propagation of coherent states can be directly deduced from the 
results about Hermite states. 

Proposition 5.10 For any zo £ 3? and any b G & > P q(J2f r ), the limit 

lim (u e (t)E(zo) , b Wick U e (t)E(z. Q )) = b(z t ) 

e— >0 \ / 

holds for any t £ R when Zt denotes the solution to the Hartree equation (130b - 

I |2 °o £— n/2 

Proof. By symmetry, one can assume m = p — q > 0. Recall that E(z ) = e~^ £ -^Zq " and start 
first with \zo\ = 1. Since t/ e (f) preserves the number, one gets 

(u e (t)E(z ),b Wick U £ (t)E(zo))= J £e- e ' —a n (e" 1 ) 
with a„ (fT 1 ) = £ m /V«(«-l)...(«-m + l)( z f"- m , C/ e (/)*fe WM C/ e (f)zf n 
By Lemma [2~4l the quantity a„ (e -1 ) satisfies 

K [e ) I < (ne) 2 |6| z^ q ar) - ( n£ ^ r\^{\i p ar,\i q ar) • 
Hence Lemma lATTI applied here with A = and \l= p reduces the problem to the proof of 

lim / a, pr , 2 i(A)— == <fs. 

The uniform estimate 
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and the pointwise convergence induced by Theorem 15 . 6 1 with z = Zo, k = [VXs + X] and e = A -1 yields 
the result. 

For a general |zo| > 0, write 

i °° (e'V n l 2 
n =o V"! 

with z' = y^y and e' = ■ By replacing the £-quantization by the e'-quantization, with 

b Wick,e' = \ ZQ \-P-1 b Wick for be & m (%-) 

H E = \z \ 2 dr E ,(-A) + \z \ 4 Q Wick ' e ' 
and (ied t u = H e u) & (ie'd t u = dF e ,(-A)u + |z | 2 Q Wick ' e 'u 

Hence the previous result applied with E'(z'q), \z! \ = 1 and the g'-quantization implies 



where z' t solves 



lim (U e (t)E(zo) , b Wick U e (t)E(z )) = \z \ p+q b(z' t ) 



^z; = -Az;+iz„i 2 (v*|z;i 2 )z; , z / !=0 =z / Z{) 



K0| 

Since this mean field equation preserves the norm \z[ \ like (l30l) does for \z t \, this implies 

z' t = \zq\~~ z t = \zt\~~ z t and \zo\ p+q b(z t ) = b{z t ) ■ 



□ 



Remark 5.11 Another proof can be obtained directly from Proposition \5.2\ after checking uniform num- 
ber estimates for U2(t,0)£l. But working in this direction is more efficient with the help of Wigner 
measures. 



6 Wigner measures: Definition and first properties 

The notion of Wigner (or semiclassical) measures is well established in the finite dimensional case. We 
refer the reader to [Bur][Gerl ][GMMP][HMR][LiPa][Tar] for details. The extension that we propose 
here to the infinite dimensional case follows a projective approach. 

6.1 Wigner measure of a normal state 

Consider the algebra of cylindrical sets ^ cy i(3f) = {X(p,E) = p^ l (E),p£ F,E £ 3§(p2f)} where 
38(p^) denotes for any pgF the set of Borel subsets of piF. A cylindrical measure /I is a mapping 
defined on 3S cy i(^) such that: 

• Ai(iT) = 1, 

• For any p G P, H P (A) = il(p~ { (A)) for A £ £%(p£if) defines a probability measure \x p on 3${p2£). 

The family of measures {jU p } p ep is often called a weak distribution. 

This notion is often introduced within the framework of real Hilbert spaces (or more generally real 
topological vector spaces). This makes no difference at this level. The real structure on iF, namely the 
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real scalar product S, is useful for the application of Bochner's theorem. For any £ G i2° the function 
z i — > g- 2m5 (z >§) i s a cylindrical measurable function and the Fourier transform of jx is well defined by 

Bochner's theorem characterizes the Fourier transform of a weak distribution. It says (see for example 
[BSZ]) that a function G is the Fourier transform of a weak distribution if and only if 

• Gis normalized: G(0) = 1, 

N 

• G is of positive type: V A, A i G(^,- — £j) > 0, 

• For any pGP, the restricted function G| p j> is continuous. 

An important point is that ^ is a separable Hilbert space. Hence the a-algebra generated by the cylin- 
drical sets, that is containing & cy i(£>f), is nothing but the Borel a-algebra, associated with the 
norm topology on 3f. A probability measure well defined on will be shortly called a probability 
measure on iF. The tightness Prokhorov's criterion (see [Sch]) has within this setting the next simple 
form. 

Lemma 6.1 (See [Sko]) A cylindrical measure \i on 3? extends to a probability measure on 2£ if and 
only if for any rj > there exists > such that 

VpeP, n({ze$f,\pz\<R n })>l-ri. 

By recalling that for any R > the ball {z G ^ : \z\ < R} is weakly compact, this can be reinterpreted by 
saying that a weak distribution jJ. extends as a Borel probability measure if and only if its outer extension 
is a Radon measure on endowed with the weak topology (see [Sch]). 

Consider a family (p e ) ee (o,e) of non negative trace class operators on such that Tr[p e ] = 1, or 
equivalently normal states & i— >■ Tr[p e <^] on the space of all bounded operators J2? (Jt?) . An additional 
number estimate assumption allows to associate with such a family, Wigner probability measures on 3f. 

Theorem 6.2 Let (p e ) ee (oe) ^ e a f am dy of normal states on J2?(J^) parametrized by s. Assume 
TrfA^p 6 ] < C§ uniformly w.r.t. e £ (0,e) for some fixed 8 > and Cg G (0,+°°). Then for every 
sequence (£ n ) n eN with lim^^o e n = the exists a subsequence (£ nk )keN an d ci Borel probability measure 
jli on 2f such that 

lim Tr[p e " k b Wey! } = lim Tr\p^- Wick ] = [ b{z) dfl(z) , 
k->°° Jg- 

for all be U peP J^ 1 {Jt h (j>T)). 

Moreover this probability measure \L satisfies \ |z| 25 d\l(z) < °°. 

Jar 

Remark 6.3 a) By introducing the reduced density matrix p e p G Jz? 1 (r s .(/?i2°)) defined for p G P as 
a partially traced operator Tr\p^A\ = Tr[p £ (A (g>/ r j p x one could consider the Husimi function 
Hp of p £ which is its finite dimensional Wick symbol. It is known that this makes a weak probability 
distribution which admits weak limits after extracting subsequences £, lk — > °°. The number estimate 
implies in finite dimension that such a limit is a probability measure. Our results say essentially two 
things: First after a proper extraction of subsequences, the family (jUp)peP makes a weak distribution, 
i.e. the convergence can hold simultaneously for all the non countable family p G P. Secondly the weak 
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distribution is a Borel probability measure. 



2.8 

b) The estimate \z\ dpL(z) < +°° will be proved in the more precise form 



l + \z\) dix{z) < liminf Tr p e »k(\+N) 



<C s <+c 



Contrary to the finite dimensional case, the first inequality is not an equality even when the right-hand 
side converges. Examples are given in Section [ 



Proof, i) The Proposition [3/7] implies 



Tr 



e ,We\'l 



Tr 



pEfoA-Wick 



for fixed b G Li pe p<^~ (^b(p^))- Hence the result is true when it is proved after considering simply 

the Anti-Wick observables. 

ii) Consider for e > the function 



e — rlsl =Tr 



G e (0=Tr[p e W(V2^ 
The positive type property and the normalization come from 
G e (0)=Tr[p e ] = l 



pS t e 2i%S{^,.)\A-Wick 



£ A ; A y G e (£-^-)=Tr 

i,j=l 



A-Wick 



>0. 



The continuity when £, is restricted to any fixed finite dimensional p^ can be written with uniform 
estimates w.r.t e G (0,£). Consider the estimate Tr [p e (l +7Y) 51 ] < C 5l with Si G (0,min(l,2<5)). Write 
for any ^,rj G 3? 



|G e (Tj)-G e (§; 



Tr 



£±&<^)- W (v/2* 5 )l£±i^ 



(/Y+l) 5 '/ 2 



(N+\) 5 '/ 2 



g 2 Ml - e 2 IS 



< 



[WiVlK^-WiVln^W+l)-^ 1 Tr[(/V+l)V 



2 o 2 o 

g 2 " _g 2 151 



We have found by Lemma [37TI two constants S\ G (0, 1) and C' s > such that 

V^tjgJT, |G e (77)-G e (^)|<C^|77-^| 5 ' [(|i?| 2 +|«| 2 ) ft/2 + l], (55) 

holds uniformly w.r.t. e G (0,e) and we recall the uniform estimate |G e (£)| < 1. Hence for any e G 
(0, e), G e is the Fourier transform of a weak distribution /i e such that 



Tr 



pEjyA-Wick 



b(z) dn e ( z ) 



holds for all b G U peP J^ 1 {Ji b (p2?)). 

iii) Actually the uniform estimate (|55T > allows to apply an Ascoli type argument after considering se- 
quence (£„)„ e N such that lim„_,.oo e„ = 0: 
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Since is separable, it admits a countable dense set j¥ = i G N}. For any £ G N the 
sequence G En {£,t) remains in {a G C,|a| < 1}. Hence by a diagonal extraction process there 
exists a subsequence (s„ k )keN sucn that for all i G N, G e (&) converges in {a G C, |a| < 1} as 
k -> oo. Set 



for all leN. 



• The uniform estimate (1551) implies that the limit G is uniformly continuous on any set .JV n 
{z G J? : |z| </?}. Hence it admits a continuous extension still denoted G in (3f,\ \^). An 
"epsilon/3"-argument shows that for any £ G lim^oo G £nj , (£) exists and equals G(£). 

• Finally G is a normalized function of positive type as a limit of such functions. 

Finally the uniform estimates |G e (£)| < 1 and |G(^)| < 1 allow to test the convergence again any 
V G ^b{p^) and to apply the Parseval identity with ft = J^~~ ! (v). From any sequence (s n ) ne ^ such 
that lim„^oc£„ = 0, one can extract a subsequence {En k )k-^<*> anc ^ ^ nc ^ a wea k distribution such that the 
limit 

b(z) dfi(z) 



lim Tr 


p Zn kb Weyl 


= lim Tr 


p £n kb A-Wick 


=/ 


tli— ><*> 








J 3? 



holds for any ft G J? (L 1 (p^ ,L p (dz,))) and therefore for any ft G •5^(5?). 

iv) The Prokhorov's criterion for jU in the form stated in Lemma 16.11 is again a consequence of the 
uniform number estimate Tr [/V 5 p e ] < C§. Fix any p G P and set d = dim/?. The operators N p = 



■■ dT{p), N p x = [l pS ® dT{I p x^)) = dT{p^) and N = dT{I) 



make a commuting family of non negative operators such that N = N p + N p ± . Thus the inequality 

de . r de 
(1 + y+A^ \'>(1 + y+ N p) 

holds for any s > 0. Hence the estimate Tr [p e A^ 5 ] < Cg implies 



Tr 



P e (l + f +N P ) S 



< Tr 



p *(l + +N ) 



<Tr 



p £ (2 + ^V) 5 <C^, 



with C' s > independent of e and p as soon as e < k- 

Let x G be a non negative function on piT, such that % = in a neighborhood of {|z| < 1}. 

For any 7? > 1 the estimates 

(1+tf 2 ^ 



(l + |z 



2^5 



X(R~ L z) < 1 



holds with uniform estimates of the left-hand side in S p ^{\ 



\dz[ 



The pseudodifferential calculus in 



with the metric 



\dz[ 

w 



, provides the inequality of bounded operators on T s (p^*) 

Weyl 



with 



(l+RyAoBnoA-Cs < 

Weyl 



< 1+C£ 



1+ z 



, B R =[ X (R- l z)] 



Weyl 



and 



with a constant C > independent of £ G (0, \) and R>1. By Proposition 
C" > independent of s G (0, j) (and 7? > 1) such that 



<C, 

there exists a constant 



i , d£ ,s 



<C'e. 
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Hence the inequality 



after tensorization with /, 



(1 +R 2 ) 8 x(R- 1 pz) Weyl < (1 +2Ce)A" 



and testing on the normal state p e yields 



(l+R 2 ) s Tr p e x(R- l pz) We - 



1 „-\Wevl 



<c'k 



with a uniform constant Cg with respect to e G (0, g) and /? > 1. After taking the limit — > °°, — ► 0, 
we get 

l {N >« } (z) rfM(z) < / = I™ Tr [p^/rVz)^! < Cg(l + /? 2 )- 5 . 

This inequality is valid for any pEP and the Prokhorov's criterion of Lemma l6.1l is satisfied. The weak 
distribution pt is a probability measure on 3f. 

v) First the function (z) is Borel measurable in 2P. Take p G P and /? > 1 and take now # € c $^\jp3f), 
such that < < 1 and #o = 1 in a neighborhood of 0. Consider the estimates 

(l +iV )5>(l +A r)5 > (i+Ar p ) 5 / 2 2o(«-V)^(l+A^) 5/2 -C p e(l+A^) 5 



> 



i + H ) Xo(ir l pz) 



Weyl 



■C p e(l+N) 1 



where the two last inequalities are again derived from the finite dimensional Weyl calculus (with a 
uniform control w.r.t. R > 1). After taking the limit — ► oo, e„ k — > 0, this implies 



J^ + I/kI 2 )**^*" 1 /*) <*m(z) 



= lim Tr 
< liminfTr 



pS 



(i + H ) zo(^Vz) 



Taking the supremum w.r.t R > 1 and then w.r.t a countable increasing sequence (p n )„ 6 N, p n G P, such 
that sup„ eN p„ = 7^, yields 

/ (l + |z| 2 ) 5 ^(z) <Cg <+oo. 
■IS 

□ 



6.2 Complex Wigner measures, pure sequences 

More general families of trace class operators can be considered by linear decomposition 

P £ = *>r+Pr+ ~ K-Pr- + &/+P/+ - »**-p/_ 
with A. e > 0, p. e > 0, Tr [p. e ] = 1 and 

A| + + A|_+A/ + + A/_<4Tr[|p e |]. 



(56) 



Proposition 6.4 Let (p e ) ee (o,e) ^ a family of trace class operators such that 



(l+A0 5/2 p e (l+A0 



5/2 



(57) 
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uniformly for some 8 > and some Cg < +°°- Then for any sequence (£ n )neN such that lim„^oo £ n 
one can extract a subsequence (Cfi^eN an dfind a (complex) Borel measure fX on 3? such that 



0, 



lim Tr[p e "kb Weyl ] = lim Tr\p e "k}^- 



Wick] 



3f 



(z) dy,(z), 



(58) 



for all b £ U peP &~ x (Jt b {p2?)). 
Moreover this measure satisfies f^-(z) d\fl\(z) < +°°- 

Proof. The decomposition (l56l > implies 



(1 +N )S/2 p e {l +N) oi, = X e +r e R+ e _ g + _ ^ ^ 



\S/2 



with 
and 



..5 



;i+a0 5/2 p. £ (i+ao 5/2 >o 



Tr 



(l+^V) 5 p. e 



Tr 



< 



(1+A0 5/ V(1+A0 



5/2 



Hence the symmetric writing with (1 + Af) 5 / 2 p e (l +N) 5 / 2 of the uniform weighted estimate ensures 
that every term p. in fulfills the assumptions of Theorem |6.2| It suffices to extract a subsequence which 
provides the convergence for all the four terms. □ 



Definition 6.5 For a family {p e ) ee ^gy satisfying 071 ), the set of Borel measures ji which satisfy (1581 ) 
is denoted ^#(p e ,£ £ (0,£)) or simply ^#(p e ). 

Such a family (p e ) e e(o,e) ( res P- a sequence (p e ") n eN) is said pure if^(p e ,£ £ (0,e)) (resp. ^(p e, \n £ 
N)j has a single element }X. 

When the family (p e ) ee (o,e) * s P ure tne limit in (I58T ) can be written with lim e ^o instead of lim, !t ^oo. 
This provides a characterization of jtf (p e ) = {/i}. For simplicity, we shall often assume that the family 
(P e )ee(o,e) * s pure, when the reduction to such a case can be done after extracting a suitable sequence. 



6.3 Countably separating sets of observables 

In order to identify a Wigner measure of \i £ (p e ) it is sufficient to test on a "dense set" of observ- 
ables. The good notion is given by the Stone-Weierstrass theorem for L 1 spaces. It can be recovered 
from the standard Stone-Weierstrass theorem for continuous functions in our case. 

Lemma 6.6 ( cf HCouV ) Let v be a Borel probability measure on a separable Banach space X and let 
{fn, n 6 N} be a countable set of bounded V -measurable functions which separates the points 

Vx,y£X,3n£N, f n (x)+f n (y). 

Then for any p £ [0,°°), the algebra generated by {f n ,n £ N} is dense in L p (X,dv). 

Since "the" Wigner measure is not known a priori, the good notion of "dense set" that we shall use is 
the following. 

Definition 6.7 A subset Q) C U pe p (^b(p3f)) is said countably separating whenever it contains 
a countable subset, *3> D S>o ~ N, which separates the point of 

Vx,y£&,3f£$> Q , f(x)^f(y). 

Proposition 6.8 Let fj,[ be a bounded Borel measure on £t? and let (p e ) e6 (o,e) be a family of operators 
which fulfills the assumptions of Definition 16.51 The two next statements are equivalent: 

1. ^(p e )={ M i}. 
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2. There exists a countably separating subset & C U pe p & 1 {^b{p3?)) such that 



Vbe@, lim Tr 



p e b Weyl 



lim Tr 



b(z) dpLi(z). 



Remark 6.9 A similar equivalence is obtained for H\ G <J?(p e ) a/ter a subsequence extraction. 

Proof. Assume G Jt{p e ). There exists a sequence {e nk )km an d a Borel measure ji. such that (158 
holds for any b G U p(E p J^ -1 ^^/?-^)- In particular this holds for any b ^ Q>: 



b{z) d\i(z) = limTr p e "*& 



We\l 



3f 



b{z) diliiz). 



The set <2l is dense in L l {3? ,d\}X\\) and in L ! (i2°, so that the above equality of the extreme sides 
extend to any bounded Borel function. This implies \x = }X\. □ 
The next examples will be useful in the application and allow to reconsider an inductive point of 
view. 

Proposition 6.10 Let (p(,)eeN be an increasing sequence of projectors in P such that sup^pg = 
and let the family of operators (p e ) e e(o,e) satisfy the assumptions of Definition 16.51 Then the identity 
^(p e ) = {jit} is equivalent to any of the next statement 

1. For all b G U^jq 5?{pi3?\ the quantity Tr[p e b Wey ' ! ] converges to f^b{z) d}l(z) as e — > 0. 

2. For all b G S^cyi{3f), the quantity Tr[p e b We> ' 1 ] converges to fg b(z) djl(z) as E — > 0. 

Proof. It suffices to notice that Uf £ N^(Pf^°)> and therefore ^cyi(^), is countably separating because 
the weak topology separates the points. □ 

6.4 Orthogonality argument 

Complex Wigner measures are especially interesting while considering the joint measure associated 
with two families of vectors (« e ) ee (o,e) an d { v£ )ee(o,e)- Introduce the notation 

P H e v=l" e )(v £ |. 



Proposition 6.11 Assume that the family of vectors (H e ) £6 (o i e) and (v e ) e e(o,e) satisfy the uniform esti- 
mates 



;i+A0 5/ V 



+ 



(l+AO 



<c 



1 



for some fixed 8 > and C > 0. Assume further that any jX G -^(pf 1( ) and any V G ^{pf n ) are mutually 
orthogonal. Then the family (pf v )e€(0,e) i s pure with 



^(pf v ,£G(0,£))={0} 



i.e. lim(w e ,& 



WeyL,e 



e^O \ 



lim ( u e , b 

£->0 \ 



,A—Wicky£ 



for any b G & (^b(p^)) and any p G 
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Proof. Assume J%{p m ) = {/i} and jK{p^) = {v} with }X L V. Take 77 > 0. There exist two bounded 
closed subset ^1 and K 2 such that 

H(Ki)>l-r] , v{K 2 )>\-T] , K l nK 2 = <d. 

Since Ki and K 2 are compact in the weak topology, Ki C CK 2 , ZK 2 open in the weak topology, there 
exists a finite covering of K\ of the form 



Kl C U {\p k (z-Z k )\<r k } 

k=l 



U {\p k (z-Z k )\<2r k }nK 2 = <d 

k=l 



with pk G P, Zk G ^ and r# > for all & G {1, . . . , A"}. By choosing for any k a function Xk £ ^oiPk^) 
such that Xk{Pk(z)) = 1 when |/7^(z — < and ^(ftz) = when \p k (z-z k )\ >2r k the sum^(z) = 
£jLj ^ Txkt(piz) d e fi nes a cylindrical function # G t 5^ cv /(^') such that £ = 1 on #1 and # = on K 2 . 
Take now any Z? G oS^ cv /(iF) and write 



V , (^)^V) I + I (u* , (b(l - x)) Wey 'v e ) 
< Ub(l-x)) w v l u e \^+\(bx) 



\Weyl v e 



From the Weyl pseudodifferential calcul we get 



2 1 u\2\^ e y^ 



{b{l-x)) Weyl u^<Tr P U{l-X?\b\ 2 ) 



where the right-hand side converges to f s \b\ 2 (l — x) 2 (z) djl{z) as e — > 0. The property x = 1 on K\ 
with (K\ ) > 1 — t] implies 



\im^ V \(b{\-x)) Weyl u £ \ 2 ^<^b\ 



and with the symmetric argument limsup g ^ | ipx) yS \ w ^ H Hit- Hence we get 



Vrj > 0, limsup 



b Weyl v e 



<2\b\ L ^ 



for any b G This implies Jt{p* v ,e G (0,e)) = {0} . 

A straightforward consequence is the next proposition. 



□ 



Proposition 6.12 Make the same assumptions as in Proposition \6.11\ with the additional condition 
^(Puu) = {/-*«} an d ^(Pvv) = {Mv}- Then the family of trace class operators (Pu+ v ,u+v)e€(0,s) sat- 
isfies 

04+ M • 



e 

u+v.u+v, 



Proof. Write simply 



u e + v e ,b Weyl (u E +V 



e ,b Weyl u e ) + (v e ,b Weyl V 



+ (u 



.£ uWeyl e 



v £ ) + (v e ,b We y l u e 



and take the limit of every term as e — > 0. 



□ 
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6.5 Wigner measure and Wick observables 



Up to some additional assumption on the state and by restricting the class of Wick observables, we 
check in this subsection that testing with Weyl, (or Anti-Wick) and Wick observables provides the same 
asymptotic information as e — > 0. 

Fix once and for all p G P, the choice of the metric g p = \dz\ 2 or g p = ^jgr- From Proposition 13 . 8 1 we 

know that the class of symbols U P £W, S £rS p %>((z) s ,g p ) and ® r ^ q ^^m,q{^') both contain all the classes 
&m,q{p^), with a good comparison of Weyl and Wick quantizations on these smaller sets. In the limit 
£ — > 0, this comparison can be carried out to any b G ^^qeN^m^i^)- 

Theorem 6.13 Assume that the family of operators (p e ) eG (o,e) satisfies 



(1+A0 5/ V(1+A0 5/2 



<c s 



uniformly w.r.t e G (0,e) for any 8 > 0. 



1. For any fixed j3 G U^p^r Spj-^z} 1 ,^), the families (P ^ p e ) E e{o,e) and (P 
satisfy the assumptions of Definition \6.5\ and 



A-Wick p£^ 



.J?{p Weyl p £ ) 



j?(p A - Wick p £ ) = {fifx , /i g (p e )} 



ee(0,e) 



(59) 



2. For any fixed j8 G ® & mqeH^ > m,q{^') the family (P Wick p e )e£(0 1 e) satisfies the assumptions of Defi- 
nition \6.5\ and 

Jt($ mck p e ) = {p> , M G ^(P e )} • (60) 
A particular case holds when the measure is tested with o = 1. 

Corollary 6.14 Assume the uniform estimate |(1 + ./V) 5 / 2 p e (l +N) 5 I 2 \^^^ < C$for all 8 > and 
further Jt(p £ ) = {/*}. 

i. 77ie equality 



lim 7r 



lim Tr 



jA-Wi'cApe 



J3(z) dii(z) 



holds when j8 G U p€ w tS&s >S pa r((zy,g p ) 
2. The limit 



lim Tr 



j6(z) dix(z) 



holds for any /J G ©:! g , eN ^~ (iT) 



Proof of Theorem 1) The relation ([27]) extends to any b G S p &((z) s ,g p ) and implies e -1 ^ 1 ^ - 
&A-Wfc*) = c ( e )Wey/ with c ( £ ) uniformly bounded in S pS -((z) s ~ 2 ,g p ). The result for p A - Wick can be 
deduced from the one for p We >' l _ 

Take p G P, 5 > (this contains the case s < 0) and j8 G S p ^-((z) s ,g p ). Let AT p = N p ^®l Y ^ p x^ and 
Af± = Ir i {pZ)®Np L &- Our assumption on (p e ) ee (o,e) and the commutations [W_j_,iVp] = [N p ± , [5 Weyl ] = 
imply for any 8 > 

(1 +A0 5 / 2 j 3 H/e - v 'p e (l +A0 5/2 = AfiA'flC 
with A = (1 +A0 5 / 2 (1 +Ag~ 5/2 (1 + )~ 5/2 
fl = (1 +N p ) S/2 p Wey '(\ +N p y S/2 - s 



A' 



x, ■ [I + N) 



R = (\+N) s+S p e (l+N) s+S and C = (1 +A0~ 5/2 ^ • 
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The factors A, A' and C are uniformly bounded operators when 8 > (and s) is fixed. The trace class 
norm of the factor R is uniformly bounded by C§ +s . Finally the Weyl pseudodifferential calculus on 
p3? implies that B = y Weyl with y(e) uniformly bounded in S p %-(l,g p ) and therefore \B 
uniformly w.r.t e G (0,e). 

Hence the family {P Weyl p 8 ) e <={o,e) satisfies the assumptions of Def.O Let ^ belong to ^(p We >' l p e ). 
After extracting the proper sequence (e n ) ne ^ such that lim„^oo£,j = 0, one can assume 



< C 



lim Tr 
and lim Tr 

n— »°° 



b Weylp 



Weyl e„ 



b(z) d/i^z) 



b Weyl £ „ 



2f 



'2f 

b(z) <ijU(z) 



for any b G y cy i(3f). But the finite dimensional pseudodifferential calculus implies b Weyl fi Wey 
(bp) We y l + 0^ {Jf) (e n ) with bp G y cy i(3f). This implies 



f b{z)d^{z) = I b(z)P(z)dn(z 

2? J St 



for all b G S^cyli According to Proposition 16. lOl fhis implies }X\ = p}X. 

2) Since the U p e¥, s eRS p3 r({z) s ,g p ) contains \J pe & \ ® a ^ qm ^m, q {p^')\ , the result is proved for any 



polynomial symbol b G £P™ q (3f) such that b = T{p)bT{p) for some finite dimensional projector p G 
Consider now a general b G ^^ q (3f) with m,q G N. By Lemma [2~4l the operator 



(1+A^) 5/2 ^(1+A^)- 



8/2-m/2-q/2 



is uniformly bounded for any 5 > 0. Since the trace class norm of (1 +_/V)^ ± ? t2 p e (l +./V)^ ± r t2 is 
uniformly bounded w.r.t e G (0,e), the family (p w,ck p e ) satisfies the assumptions of Definition 16.51 
Introduce now an increasing sequence (pe)i e ® of F such that sup fcN pe = I and consider for I G N 

Pt(z)=P(p t z) , Pt=pf q o~bopf>». 

Since /3 is a compact operator, the finite rank operator pg converges to p in the norm topology in 
JS?(V m \J q 2?). The uniform estimates 



(P-p e ) Wick (i+Ny m/z - q/2 

? \ m/2+q/2 



<c 



jSf (V w iT, V" JQ 



and the convergence 

VZ?GXy/(^), limTr 



(|j3(z)| + |fc(z)|)<C with limft(«) = P(z), 



foWeyl pWick p£„ 



b(z)Pi(z) dn{z) 



after extracting a sequence {e n ) ne m, lim n ^oo£„ = 0, with J^ (l + \z\ 2 )' n / 2+q l 2 dfx{z) < +°°, lead to 



ybe^cyii^), limTr b Weyl p Wick p e " 



(z)P(z) d^{z). 



□ 

The previous results provide the behaviour of lim e ^ Tr [p Wick p e ] for J3 G ©^ GN ^m, ? (^) when 



^(p e ) = {jj,}. The next result checks the other way. 
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Proposition 6.15 Assume that the family (p e ) ee (o,e) satisfies rt57D and that for any C > there exist 
Kc > smc/i that 

£M! z ' rlivVl - Kc< ~ 

/jo/<is uniformly w.r.t £ G (0,e). Assume that there exists a Borel measure jj, such that 



lim Tr 

e^O 



.2 



6(z) d\x(z) 



holds for any b G (Bm S q&m,q{^')- This implies 

Jt{p e ) = {ll}. 
Proof. It is enough to prove the following statement: 



limTr[W(£)p f 

£->0 



e^ s ^d]X. 



It is done when the right-hand side of 
Tr[W«)p«] = 



£J^i!! Tr 



n=Q 



2"n\ 



Wick 



\ 1^1 



(61) 



is proved to be an absolutely convergent series, uniformly w.r.t. e G (0,e). With 
Tr[W(£)p e ] = limTr[l^)l [0 , M] (/V)p e ] 



/>/25(g,z)' 



Wick 



1[0,M1(^)P £ 



(62) 



and 



Tr 



Wick 



<M n 



Wick 



-n/2 



Iv^l 



with M„ = Tr[(l +A^)"p e ], Lemma EH implies 



Writ 



-n/2 



(l + 2ygg+7jef n! 
" ^ P N (te + l)»/2 (j - £ + l)»/2 [„/ 2 ]! 



i! 



< 



[n/2]f 



This leads to 



7J=0 



2"n! 



Tr[/j„ 



{ 1^1 



Writ 



* I n^rr M " <o ° (63) 



uniformly w.r.t. e G (0,e) and M > 0. Hence we can take the limit M — > oo inside in all the terms of 
(l62l . This leads to (I6TI ) with a uniformly absolutely convergent series in the right-hand side according 
to (I63T ) and our initial assumption. 
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Thus the sum and the limit as e 



can be interchanged in (I6TT ): 




) 



Wick 



limTr[W(£)p e ] 




The last equality follows owing to the dominated convergence theorem and 




d^ = Hmy-Tr[p £ dr(\ pS ) k ] 

e ^ u k=o K - 



CO 



for any 8 > and any pgP. This completes the proof. 



□ 



7 Examples and applications of Wigner measures 

7.1 Finite dimensional cases 

The first examples are given by Theorem l4.2l 

1. For any z G 3f the family of operators p e = \E(z))(E(z)\ has a unique Wigner measure 



2. For any z G ^ and any m G the family of operators p e = |z®^ E "^(z®^! with |z| = 1 and 
lim g ^o = 1 nas a unique Wigner measure 



3. In case 1) and 2) the convergence can be tested with Weyl, Anti-Wick of Wick observables ac- 
cording to Proposition 16.41 and Theorem l6.13l 

Beside the explicit calculation of Theorem 14.21 these results can be considered through an inductive 
approach since E(z) or z® n lie in Y s (Cz). The natural extension comes from Proposition 16. 1 01 - 1 ) with a 
proper choice of the first term in the increasing sequence (pe)ieN- 

Proposition 7.1 Assume that the family (p e ) ee (o,e) satisfies the assumptions of Definition 16.51 Assume 
further that there exists a finite dimensional space poGP such that 



for all e G (0,e) with p pa G Jz? (T s (po^f)). Then the Wigner measures 0/(p e ) ee (o,e) are given by 



j$r(\E(z))(E(z)\, ee(0,e))={8 z }. 




P e = r(p ) P r( Po ) = P e po ®\a)(a 
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7.2 Superpositions 

Two kinds of superpositions can be considered : 1) convex or linear combination of trace class operators; 
2) convex or linear combination of wave functions. The first one is the simplest. 

Proposition 7.2 1. Let (M, n) be a probability space. Les (p e (m)) ee( - ^ meM be a family of opera- 
tors such that 



l+N) S/2 p e (m)(\+N) 



8/2 



<C s (m) 



for %-almost every m £ M with Cg G L 1 (M, dn) for some 8 > 0. Assume further .<#(p e (m) , £ G 
(0,£)) = {jU(m)} for %-almost every m€M, then the family (JMP E ( m ) < ^ ?r ( m ))e66(0,e) satisfies 
the assumptions of Definition \6.5\ and 



J p e (m)dn(m), £6(0,£)^=jy /i(m) cf7r(ra) j . 



2. Any bounded Borel measure on 3f can be achieved as a Wigner measure. 
Proof. 1) Set p e = f M p e (m) d%{m) and write 

(1+A0 5 V(1+A0 5/2 „ < / C s (m)dn(m). 
Then apply Lebesgue's convergence theorem to 



Tr 



Tr 



M 



b Weyl p e (m) dn{m). 



2) After reducing the problem to the case when /J. is a Borel probability measure on apply 1) with 
M = £¥,7i = lJ,,m = z and p e (z) = |£(z))(£(z)|. □ 
The second type of superposition requires an orthogonality property. It is given by Proposition ^. 121 
Here are a few examples 

1. Take u\ = E(zg) for £ = 1, . . . ,L, with L G N fixed, and set u e = L _1 / 2 ^f =1 u\ . When the zi are 
distinct, the family ( |w e ) {w e |) £e (o e) has a unique Wigner measure 

^{\u £ )^\) = h-^8 z \. 



2. Take for any I G {1,...,L}, u\ =zf ke with |^| = 1 and lim£A: e = 1. The family (\u e }( 



has a unique Wigner measure: 

Jt{\u e ){u e \) 

.e _E{z)+\z®'- 



L f 27C ) 
(2^L)- 1 £y () 8 e ,e 7( dd\ 



3. For z G i2° and w e = — - with \z\ = 1 and lim £ ^,QEk E = 1> the family (|« e )(w e |) eG (oe) has a 
unique Wigner measure: 



^(| M e )( M e |) 



1 /' 2?r 



4. All this examples can be tested with Weyl, Anti-Wick or Wick observables according to Proposi- 
tion [63] and Theorem 16. 131 
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7.3 Propagation of chaos and propagation of (squeezed) coherent states 

Let us go back to the example of Section [5] where U e (t) = e^ He with H e = dT(-A) + Q Wick , Q = 
\V{x\ — x 2 ) and z t solution to id t z t = — Az + (V * \z.t\ 2 )z.t Theorem 15.61 Proposition 15. 101 and Proposi- 
tion [6j3] imply: 

1. For any z G 2t with \zo\ = 1, the family {\U £ {t)zf e ) {U e {t)zf e |) ee m with lim e ^ £^ e = 1 is 



pure with 



ur (\u £ (t)zt)(Ue(t)z^\) = jfV* de}=^ (|zf E )(zf E |) 



2. For any zo G the family (|J7 e (f)£'(zo))(t/e(0 £ '(^o)|)ee(o,e) is P ure witn 

^(|J/ e (0£(zo))(f/ e (^(^o)|) = {5 Zf } = ^(|£fe))^fe)l)- 

These results are derived from the results for product states after testing with Wick observable (any 
b G ®n^q^ > m ,q{3f)) ■ Actually it is possible to recover the second one directly from the Hepp method. 
For any b G S"cyi(3f), Proposition 15.21 implies 



lim Tr 



b w * (\U e (t)E(z ))(U e (t)E(z )\ - \wC^z t )U 2 {t,0)Q.)^{—Zt)U 2 {tm\ 

\ IS IS 



0. 



By the finite dimensional Weyl quantization, the second term equals 

(u 2 (t,o)nM--zt) We) 'u 2 (t,o)n). 

And it suffices to check that the family (\U 2 (t, 0)Q.) (U 2 (t,0)Q.\) ee ( 0e - ) admits the unique Wigner mea- 
sure do- This is a consequence of Lemma [531 which first says \N k U 2 (t,0)Cl\ ^ < C& for any k > and 
then lime-tf (02(/,O)Q, b Wkk U 2 (t,0)Q.) = when 6(0) =0. 



7.4 Dimensional defect of compactness 

In the last example the mean field propagation of Wigner measure attached with U e (t)E(zo) can be 
proved directly without using the result on Wick observables. As a corollary, this provides the result for 
Wick observables b Wwk when b G according to Theorem [6-131 The result for a general 

b G (Bm g q&m,q(^) is still true but comes from a direct proof or from Proposition 15. 101 
A natural question is whether the result of Theorem 16. 131 can be extended to any observable b w with 
b G (Bm^^m^i^)- The answer is no, because in the infinite dimensional case there can be some defect 
of compactness w.r.t to the dimension variable. 

Here is a typical example. Consider a family (z e )ee(0,e) such that z e converges weakly to 0. There 
exists a constant C > such that |z e | < C for all e G (0,e) and the family (E (ze)) e€(0,e) satisfies the 
assumptions of Proposition 16.151 The Wigner measures }i G ^{\E{z e )(E{z e )\)) are determined by 
testing on any b G &™ q {3f). But Theorem [42] says 

(E( Ze ),b Wkk E(z £ )) = b(z £ ) = (zPMT) ■ 

When m + q > 1 the operator b is compact, the right-hand side converges to as e — > 0. According to 
Proposition 16. 151 this implies 

^(\E( Ze ))(E(z E )\) = {So}. 
Meanwhile testing with N = dT(I) = (\z\ 2 ) W ' C implies 

(E(z e ),NE( Ze )) = \z e \ 2 
where the right-hand side can reach any possible limit in [0,C]. 
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7.5 Bose-Einstein condensates 



The thermodynamic limit of the ideal Bose Gas presented within a local algebra presentation in [BrRo] 
can be reconsidered by introducing a small parameter e — » 0. Namely, the large domain limit where 
bosonic particles are moving freely in a domain A, with volume |A| — ► °°, can be formulated with 
|A| = i and e — > 0. For a fixed particle density the total number of particle is 0(j) coherent with a 
mean field approach. Before considering any dynamical problem, Wigner measures of e-dependent 
Gibbs states bring some interesting presentation of the Bose-Einstein condensation. 

Consider the Laplace operator H = —A x on the £-dependent torus M. d / '(e~ x l d 7L) d with spectrum 
a (H ) = {s 2/d \27tn\,n G Z d }. The one particle space is jT £ = L 2 (R d / '(£- l / d Z) d ) and the bosonic 
Fock space is M' e = T s (3f e ). For the inverse temperature j3 = ^ > and a chemical potential }X, the 
Gibbs grand canonical equilibrium state is associated with the operator e~/^ r (#o-M/) = p^ e -^(i?o-^)) s 
which is trace class if and only if pL < (see [BrRo, Proposition 5.2.27]). This Gibbs state on Jz? (J^ e ) 
is given by 



(O e (A) =Tr[p e A] with p e 



1 



Tr [r(e-P(«»-M))] 



T(e 



At <0. 



It is convenient to introduce the parameter z = e?^ and this Gibbs state restricted to the CCR-algebra 
(the C* -algebra generated by the Weyl operators W\ (f), f & 3f £ ) is the gauge-invariant quasi-free state 
given by the two-point function: co e (a\(f)a\(g)) = (g ,ze~P H °(l — ze~P H °)~ l f) . The index i means 
that the CCR are written at this level in their initial form: [ai(g),a\(f)] = (g, /). This is proved in 
[BrRo, Proposition 5.2.28] with the straightforward rewritting 

co e (W,(f)) = exp [-(f, (1 +ze-^)(\-ze-^Y x f)/A 
The mean field analysis consists here in introducing a(f) = e l l 2 a\ (f) and W(f) = Wi (e l l 2 f): 

co e (a*(f)a(g)) = s(g,ze-P H Hl-z.e-^)- 1 f) 

co e (W(f)) = exp [s(f, (1 +ze-^)(\-ze-^)- l f)/A . 

Further a rescaling motivated by the observation of the phenomena on a large scale, is implemented 
with f(x) = E l l 2 (p(e l l d x) = D e (p. After conjugating with the unitary transform T(D e ) : Jif = T s (^) — > 
Jff' e = T s (3f e ), with 3f = L 2 (M. d /Z d ) we are led to consider the asymptotic behaviour as e — > of the 
normal state 

p e = T(D e )*p e T(D e ) = - p- „* ,^ T(e-^ Vd ^) 



which satisfies 



Tr [T(e-^- £2/dA -^)] 



Tr[p £ W(f)} = exp -^/,(l+z/ e2/rfA )(l-z^'" A )-7).r 



j3e 2 / d A\-l 



_e if a 

? 4 u 12 exp 



^(f,ze^ /dA (l-ze^)- l f) s 



Tr[p £ a*(f)a(g)] = e(g,ze^'\\-ze^)- 1 f) s . 

The above expressions are explicit after the decomposition in the Fourier basis / = ^nez d fn e 
any element / e Jf . For a given z < 1 and j3 > the rescaled particle density is given by 



liitn.z 



Of 



ez 



z.e 



-^e 2 l d \2nn\ 2 



ez 



+ e V 

1-7 V (1 - 7 p-Pe 2/d \27m\ 2 \ 1-7 



+ V e (j8, Z ) 



(64) 
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One checks easily for e' > e and z' < z < 1 



MjM<v e (j3,z) e ^°vo(JM = 
and Vee[0,l), v e (j3,z)>v e (/3,z') 



ze 



-j}\2mi\ : 



1 — ^-^|27t»| ; 



du 



Here comes the discussion about the Bose-Einstein condensation. In dimension d > 3 (this restriction 
may change with an alternative Hamiltonian Hq = X{D X )), the quantity 



voOm; 



j3|27tw| : 



-P\2ku\- 



du < +oo . 



is well defined. 

We focus on the case d > 3. 

The previous discussion imply 



V £ >0,Vze(0,l), v e (j8,z)< v (j3,l) 



while any total density can be achieved by (164b . The Bose-Einstein condensation occurs while consid- 
ering the limit e — > with the constraint -j^f- + v e (ze , /3 ) = V with /3 > and v > fixed. There are 
two possible cases: 

•v < Vo(j8, 1): Then lim e ^ z e = z < 1 and lim e ^ j^r = . 

•V > Vo(j8, 1): The inequality v - Vo(j3, 1) < < V leads to z e = 1 - — ^ ^ +o(e) . The propor- 
tion 1 — Vo(j8, l)/v of the gas lies in the ground state n = of the one-body Hamiltonian. This is 
the Bose-Einstein condensation phenomenon. 

It is interesting to reconsider this limit e — » with j8 > and v > fixed (d > 3) within the Wigner 
measure point of view. This is possible owing to the explicit formula 

-PE 2 ^\2Kn\ 2 



Tr 



p e W{Vlnf) 



,-eAf\\ 



exp 



-en 



I \fn\ 



z e e 



(\- Ze e-P eVd \ 27t "\ 2 ) 



(65) 



where / = Y<neZ d fn^ nx - Remember that the charactistic function of Wigner measures are determined 
after considering the limit e — > of the above expression for any fixed / S 3f. Hence the problem is 
reduced to the application of Lebesgue's theorem in the argument of the exponential. 



For any n ^ the quantity 
we get 



converges to as e — > because d/2 < 1 and z e < 1. Hence 



limTr 



p £ W(Vlnf) 



lim exp 

e^O 



enhe . 2 
T^ l/o1 



With the constraint j^- < v < +°°, there are two possibilities 

• First lim e ^ = implies v < V (j8, 1) and J?(p e ) = {5q}. 



The second case lim e ^o = V — Vq (j3 , 1 ) > implies 



limTr 



p e W(V2nf) 



,-^(v-v (j3,l))|/o|- 



_^(v-vo(j3,l))|</,l>| 



Hence the Wigner measure of the family (p e ) e >o equals y v x 5o on 3? = CI x {1}^ where y v is 
the gaussian measure 



7v(zi) 



(7r(v-v (j8,l)) rf / 2 



Zi EC. 
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Our scaled observables can measure asymptotically only the Bose-Einstein phase in a non trivial way. 
The rest of the state provides the factor 5o- While testing with the observable (Izl 2 )^"^ = N, the dimen- 
sional defect of compactness phenomenon already illustrated in Subsection !7.4l occurs again: only the 
density of the condensate remains. 

Remark 7.3 i) It is possible to consider various dispersion relations Hq = X(D X ) and the discussion 
about the dimension may change. Other boundary conditions (here periodic boundary conditions 
are considered) and the discussion about the convergence oflim e ^oz e = 1 may change a little 
bit. We refer the reader to [BrRo] for the case of Dirichlet boundary conditions. 

ii) From (|65T > it is possible to consider the limit for any fixed f G 3f as £ — > with various behaviours 
of Ze- This provides asymptotically a weak distribution. But the uniform tightness assumption 
Tr [p e (l +N) S ] < C is not satisfied. The scaling has to be adapted differently to the dimension 
d = 2 or d = 1 by taking care of the singularity at the momentum 0, in order to allow a non trivial 
Wigner measure in the thermodynamic and mean field limit. 

7.6 Application 1: From the propagation of coherent states to the propagation of chaos 
via Wigner measures 

In the previous sections we showed how the propagation of (squeezed) coherent states can be derived 
from the propagation of Hermite states or directly via the Hepp method. The Hepp method is very 
flexible (see [GiVe] for example) and therefore it is interesting to know whether a result for coherent 
states provides an information for product states or more general states. Here is a simple and abstract 
result which relies on some gauge invariance argument. 

Theorem 7.4 Let U e be a unitary operator on .Jtf possibly depending on £ G (0,£) which commutes 
with the number operator [N,U e ] = 0. Assume that for a given z G ^ such that \z\ = 1, there exists 
Zu G 3? such that 

^(\U e E(z))(U e E(z)\) = {8 Zu }. 

Then for any non negative function <p G L l (R,ds) such that f^(p(s)(\ + \s\) s ds < °°for some 8 > 
and J R <jp(s) ds = 1, the state 

oo 

pl = Y,z y ' 2( P(z il2 (n-e- x ))\U e z m )(U e z m \ 
satisfies the conditions of Definition \6.5\ and 




Proof. Owing to the relation 

r(e^ i9 )b Weyl r(e ie ) = e - idN b We y' e WN = b(e~ ie .) Weyl . 

Our assumptions imply 

JK (r(e ie )\U e E(z))(U e E(z)\r(e- w )) = 8 e ,e Zu 

for any 6 G R. The assumptions of Definition 16.51 are satisfied because U e preserves the number. After 
taking the average w.r.t 6 G [0,27r]: 

1 r 2n 

° £ = 2lil r (e W )\U e E(z))(U e E(z)\r(e- ie ) dd 
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this implies 

where the right-side is an extremal point of the convex set of Borel probability measure which are 
invariant after the natural action of S l on J?: S l x Jf B (y,z) — > Jz G -2?°. 
Again the commutation [t/ e ,iV] = and the expression (@]) for £J(z) imply 

a e = (2k)- 1 U E \T(e w )E(z))(T(e w )E(z)\U*de 
Jo 

fin 

= (2k)-' / U e \E(e w z))(E(e w z)\U* e dd 
Jo 



oo -I 

e e 



n=0 

For any ft G o^ C3 ,;(5°), the quantity 



£ -^-([/ e z® n , b Weyl U e z®) = Tr 



n=0 



e"n! 



converges as £ — ► to (27T)- 1 b(e w z.u) dd . By LemmalA~T1this implies 

2 

Vbey cvl (&), lim / flfe-v^-nCe- 1 ) -= = (2n)- 1 / d0, 

e^OJu 1 J V27T v/o 

where [t] is the integer part of t G R and 

a n (e- l ) = (U e z rm ,b We y l U e z® n ). 
Call 7 the Gaussian measure e 2 -24= on R. For any finite subdivision J 1 = \h 1l\ of R = h U . . . U 



/27r 

with intervals, the states 



satisfy the assumptions of Definition 16 . 5 1 with the gauge invariance 

T{e id )af l T{e- ie ) = af t . 

Moreover the state 

<f= [ \u £ z^ e ~ 1/2s+e ^){U eZ ^ e ~ y2s+e ^\ dy(s) = £ y(k)a 1 



is a finite barycenter of the af with a unique Wigner measure (2k) f 8 e ie Z(J dd. Since J" is finite (or 
countable), from any sequence (of") with lim„^oo£„ = 0, one can extract a subsequence (e nk )teN such 
that 

^(o*" k ,keN) = {v ( } . 



Since the measure /ly is an extremal point in the convex set of gauge invariant probability measures, all 
the Vi have to be identical to Since this holds for any sequence (£ n ) ne N> we nave proved for any 
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interval / = (a,j8) with a < p\ ^#(o/,e G (0,e)) = {/if/}. 
Now take i/a G L 1 (R, 7) and consider the state 



2y 



dr( J ) = £r(/ i )<r|. 



If there exists 5 > such that / K (l + |^| ) 5 v(i') dy{s) < +°°, the family (o_V)ee(o,e) satisfy the assump- 
tion of Definition 16.51 Let (£ n )neN be a sequence such that J((d$,n G N) = {v}. Fix b G 
The function y/ can be approximated in L l (W,dy) by i//" c G ^ C °(R). After choosing a finite subdivision 
such that the diameter of any Ig intersecting the support of \j/ c is bounded by A one gets 



Tr 



b Weyl q^ 



Tr 



7(1 



f 



a(Yc)^+\\Y- WcWv 



where (0{\\r c ) is the continuity modulus of \j/ c . Hence the right-hand side can be made arbitrarily small, 
uniformly with respect to £„, while we know that the second term of the left-hand side converges when 
\jf c and are fixed. We have proved 



b{z) dv(z) 



lim Tr 

n— 



jyWeylpE,, 



b(z) dixuiz) 



for any b G y cy i(3f) and this proves v = /ij/. Since this holds for any v G ^(cV), we obtain 
The result for p£ comes from 



< 



<P-IX 1/2 ( / 9(0 *) 1/| 



l,ds) 



with = [e l l 2 k - £- l l 2 ,e l l 2 (k+l) - fT 1 / 2 ] and = 9(j)v^we T . The condition / K (l + |s|) 5 
(p(s)ds < +00 ensures that ^#(p|) is well defined. □ 



7.7 Application 2: Propagation of correlated states 

This a simple application of the orthogonality of Wigner measures combined with the results of Subsec- 
tion [73] 

Let H e = dF(—A) + Q Wwk be the Hamiltonian studied in Section [5] and let z t denote the solution to 

id t Zt = — Azr + iy * |zr| 2 )zr- The family of integers (& e ) es (o,e) * s assumed to satisfy lim^o^e = L 

1. Let zo.e G 3T, I = 1, . . . ,L, satisfy |zo/| = 1 and set u e = L -1 / 2 J^ =1 zg e , w e (» = e'^Ue. At 
any time ? G R the identity 

jZ(\u e (t))(u e (t)\) = I (27iL)- lJ £J o 8 e ,e Zt( dd 

as soon as ziji ■ ■ ■ ,Z£,t are linearly independent. In particular this holds for any t G R when L = 2 
and zo.i and zo,2 are linearly independent. 

2. Let zo G 2 satisfy |z | = 1 and set u e = 2- x l 2 z^ +2- i / 2 E(zo) and u e (t) = e'^Ue. Then 

^(\u £ (t))(u e (t)\) = S^8 Zt + ^fJ 8 e , Ht dd } . 

3. Moreover the convergence can be tested with Weyl, Anti-Wick and Wick operators according to 
Theorem l6.2l and Theorem l6.13l . 
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A Normal approximation 

We prove a technical lemma which is a slight adaptation of the normal approximation to the Poisson 
distribution. Recall that for all — oo < « < /3 < oo we have the well known fact: 

2 

A" fP e~^r 

lim Y — e - x = -=ds. (66) 

Lemma A.l Let {a„(X)} ne x,x>o be a family of complex numbers with a n (A) = ifn < 0. Assume that 
there exist }1 6 N and > such that: 

I n \ -P 

sup \a n (A)\( Y ) < c n- 
«6N,A>0 X/L/ 

r/je« ?/ze equality 

2 

°° A" Z 1 e~ ^ 

a£^ e_A = ni^i^ (67) 

/ioWs w/zeTiever owe of f/ze fwo limits exists. 



Proof. Notice that both the series and the integral in (1671) are absolutely convergent for finite values of 

i n_\ 



A. By hypothesis a„(A) = a n (X)(j) ^ are bounded and moreover they satisfy 



} im _ E TJ^W (l - = 0, (68) 



since we may bound uniformly for A large each of the terms inside the sum and the integral respectively 
by 

CiV— e~ A n' 1 <c2, and / Lr^ ^=Ldj < Cl VA > 1. 

Therefore there is no restriction if we assume all a„ (A ) bounded by 1 since if we prove (l67l) for a„ (A ) 
then it holds for a„(A) by the limits (l68l) 



For all /i > there exists a < j8 such that 

2 2 

/ —=ds<hl, / —=ds<hl. 

Jp V27T •/-«> 

Now by d66l ) we have 



hm > — e = / —=ds, hm > — e = / —=ds 

p „ «! 7j3 V27T l-"",^ a «! J-ooy/lK 

I" vX 

Therefore there exists Ai such that for all A > Ai we have 
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Let denote I a ,p(X) = J a^ s+ ^(X)-j=ds. We obtain for all A > Ai: 



2 



1^- ^w-y^vr^w^* 



< 



£ ^-^ n (A)-/^(A) 



+2/i/3 (70) 



Using the Stirling formula there exists X2 such that for all X > X2 we have 



^ n! V27tn(n/e) n 



)] a n (A) 



< A/9- 



This yields the following estimate 



I 

n—X 
s/I 



1 ^<P(f)_ e -(^) 2 /2! 



+ 



£ — i^^a n {X)-I a p{X) 



+V12, (71) 



L a,pW 



where <p(x) =x — 1 — xln(x). To complete the proof one needs to estimate infinitesimally the two terms 
in the r.h.s. of the above inequality. Notice that by means of Riemann sums we have 



1- I 



(^) 2 /2 fP e - s y 2 



e 'vx 



VlnX 



1 

J a 



2% 



ds. 



(72) 



We have 



J_^<P(f)_^) 2 /2 1 = 



\j2%n 



e -<w) 2 / 2 



„_i „ \!2%n 

a<s 7T<P 



1], 



where <p(x) = x— 1— xln(x) + (x— l) 2 /2 which is an increasing function null at 1. Therefore one 
obtains 



' fg A9(i)_ g -(^) 2 /2i 



„_i „ y/2nn 



</'^* [^<A +1 >-1], 
./a V27T 



(73) 



with a r.h.s. converging to when A — > 00 since lim^^ e' 1 ^^ vx +1 ^ = 1 , which we bound by A/12 for A 
larger than a given A3. One can obtain the estimate 



e -(^) 2 / 2 



using the fact that 



I 

n—X 



e -(^) 2 / 2 



+ h/\S, 



a<n _x K p V2la 



(S=k)J- + 1 



< A/18, 



(i) 
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since lim^_ > „(l) = and the sum is uniformly bounded by (Equ.l72l). By splitting the integral in I a p (X) 



over the intervals [^jp ) one can snow 



that 



-A 

a 



2% 



ds 



< A/18. 



This yields 



L M (A)<V9+ I [ 



(Sr) 2 /2 r "+l-X _,2 /2 



n-X 

a 



ds] 



(74) 



with a r.h.s. converging to when X —>■ °° which we bound by h/lS for X larger than A4. Combining 
the estimates (T7Tb . (|73T > and (1741 with d70l ) we obtain that for all h > 0, there exists Ao such that for all 
A > Xq we have 



00 y n 

L^n(X) 



2n 



ds 



<h. 



This gives the claimed result. 



□ 
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